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ABSTRACT 


The stability of a dual-spin satellite system during the momentum 

wheel spin-up manuever is treated both analytically and numerically. 

The dual-spin system consists of: a slowly rotating or despun main- 

body; a momentum wheel (or rotor) which is accelerated by a torque 
motor to change its initial angular velocity relative to the main part 
to some high terminal value; and a nutation damper. A closed form 
solution for the case of a symmetrical satellite indicates that when 
the nutation damper is physically constrained from movement (i.e. by 
use of a mechanical clamp) the magnitude of the vector sum of the 
transverse angular velocity components remains bounded during the 
wheel spin-up under the influence of a constant motor torque. The 
analysis is extended to consider such effects as : the motion of the 

nutation damper during spin-up; a non-uniform motor torque; and the 
effect of a non- symmetrical mass distribution in the main spacecraft 
and the rotor. An approximate analytical solution using perturbation 
techniques is developed for the case of a slightly asymmetric main 
spacecraft. From the numerical results for the case of small mass 
asymmetry the system behaves similarly to the case of a symmetrical 
satellite; whereas for large asymmetry one component of the transverse 
angular velocity has an amplitude much greater than the initial value. 
For the case of an asymmetrical, spacecraft when the nutation damping- 



/ 


is activated during spin- up, a decay of the amplitude of the transverse 
angular velocity vector is noted. YThen the effect of the misalignment 
of the main spacecraft (spin) principal axis from the geometrical 
(polar) axis of symmetry is considered, a problem of stability could 
arise due to the large initial amplification of the system nutation angle. 


11 



NOMENCLATURE 


A,B,C 


A' ,B' ,C' 


A,B,C 


C(x,b) 


f 

F c 

F cor 

1 ( 1 ) , 1 ( 2 ) , 

X(3) 


lR i 



moments of inertia about X,Y,Z axes, respectively, for the 
main body 

composite moments of inertia about X,Y,Z axes, respectively 
including main body and rotor 

composite moments of inertia about X,Y,Z axes, respectively 
including' main body, rotor, and damper 

coefficients occuring in the solution of the differential 
equations 

Boehmer integral appearing in approximate perturbation 
solution 

spacecraft center of mass offset 
centrifugal force 
Coriolis force 

functions appearing in the particular part of the approxi- 
mate perturbation solution 

moment of inertia tensor of satellite main body i,j = x,y,z 
moment of inertia of rotor about X,Y,*Z axes, respectively 
(i = x,y,z) 

moment of' inertia of the pendulous damper(s) about X,Y,Z 
axes , respectively 
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K 


k 

a 

L. 



M 




P 

Q 



R. 

1 



= restoring spring constant of the torsion wire support 
= constants appearing in solution for the case of a symmetrical 
spacecraft and rotor without damping 
= damping rate constant 
= height of damper plane above X,Z plane 

= applied external torque about the quasi- coordinate axes of 
symmetry 

= the motor torque , 

= rotor torque 

= the mass of the main satellite 

n 

= the total system mass = M + E m. 

i i 

= the pendulum end mass 

= the inertia (applied) torque action on the rotor 

= torque due to centrifugal force about the damper hinge 
point 

= torque due to coriolis force about the damper hinge point 
= U'-B') Cu> (0) B* + I R s(0)] /B’ 

y 

= (B'-A'J I R /B' - I R 

y y 

= the generalized forces occuring in the <j>^ equation 
= the reaction of rotor on main body 

= the reaction of main body on rotor 
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r l 

s (x,b) 




s 



T 

t 


V M/ cm 

■ -V^ia’j/ cm 

V 

X,Y,Z 

x 

a 

e 

rCb.ix) 

Y 

♦i- 


= the distance from the nominal spin (Y) axis to the pendulum 
hinge point ^ 

= the length of the pendulum 

= Boehmer integral appearing in approximate perturbation 
solution 

= spin rate of rotor relative to main body 
= kinetic energ y of main body and damper 
= kinetic energy of rotor 
= kinetic energy of the system 
= time 

= velocity of main part relative to the center of mass 
= velocity of ith mass relative to the center of mass 
= potential energy associated with restoring torque effect 
= principal axes of main satellite 
= |Qc/A' | (t+a) 2 /2 

“ I xy = -“yz 
= Sz 

= incomplete Gamma function with complex argument 
= 

= damper displacement angle 

® angular velocities about the X,Y,Z, axes, respectively 
Ci = x,y ,z) 
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INTRODUCTION 


Haseltine 1 , Likins 2 , and Bainum et al 3 have investigated the motion 
of spinning^ satellites with nutation damping together with attitude 
stability criteria. A dual-spin spacecraft may be considered to 
consist of two parts constrained so that the relative motion between 
them is restricted to a common direction fixed in both bodies. Such a 
system can resist the effects of external torques because of the com- 
bined resultant momentum of the system, even though one of the parts 
may be rotating very slowly (or even with zero inertial angular 
velocity) . A basic result of dual-spin attitude stabilization studies 
is that the inertial attitude of the spin axis of a freely spinning 

W'' 

passive dual-spin vehicle may be stable, in the presence of a properly 

positioned internal damper, regardless of the mass distribution of the 

spacecraft. This means that a dual-spin system, depending on the spin 

rate and the amount and location of the dampers, may be stable in spin 

3 

about an axis of minimum moment of inertia. 

Since 1962, when the feasibility of a dual-spin satellite system 
was demonstrated with the Orbiting Solar Observatory I (OSO-I) , the 
dual-spin concept has been applied to other satellite systems including 
the Small Astronomy Satellite - A (SAS-A) , the TIROS-M Meterological 
Satellite, the TACSAT Satellite as well as the advanced versions of the 
OSO-series. The stability theory and design of dual-spin satellites 
with various types of nutation damping systems was described at a 
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symposium on the attitude stabilization and control of dual-spin 
satellites'* . Considerable attention has also been made in the recent 
literature to the dynamics and stability of various types of dual-spin 
systems with and without energy dissipation on both the high spinning 
part (rotor or momentum wheel) and the despun portion (or main body) . 
The application to dual-spin systems of heuristic or "energy sink" 
arguments has indicated that, even in the presence of energy 
dissipation, spacecraft of this type can be stable in spin about a 

centroidal principal axis of either maximum or minimum moment of 

. . ? 
inertia. 

In the paper of Mingori^, an alternate procedure for analyzing 
the motion of dual-spin satellites is presented. This procedure, which 

* ■ •• r'-j • •• ' 

involves Floquet theory, furnishes a means of obtaining precise and 
accurate stability information for a broad class of dual-spin systems, 
provided that energy dissipating mechanisms and internal moving parts 
are specifically included in the mathematical model. 

In all of these previous analyses of dual-spin systems it was 
assumed that the rotor spins at a constant relative angular rate with 
respect to the main part. Of interest in this investigation is to in- 
clude the effects of a variable rate of relative rotation such as may 
be encountered during the deliberate spin-up maneuver of the momentum 
wheel. This can be accomplished by an on-board torque motor which 
accelerates the wheel until the desired relative spin rate is obtained, 
at .the same time the main part is decelerated as the momentum is trans- 
ferred between the main part and the wheel. 
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A dual-spin Small Astronomy Satellite (SAS-A) was designed and 
developed for NASA Goddard Space Flight Center by the Applied Physics 
Laboratory and was launched in December 1970. The satellite was 
designed to scan the entire celestial sphere to determine the location 
of X-ray emitting sources relative to the fixed position of the stars. 
No serious attitude stability problems were encountered during the 
relatively short time required for wheel spin-up. Because of the 
experience already gained with this operational system, it was 
selected as a representative dual-spin system model for the present 


analysis. 



II. ANALYSIS 


The elements of the SAS-A attitude control system are shown in 
Fig. la. The satellite is comprised of three parts: (1) the primary 

part of the satellite, assumed to be essentially a right circular 
cylinder where the nominal spin axis is the body Y axis, (2) the 
smaller rotor or momentum wheel which is assumed to be connected near 
the center of mass of the primary part, and whose spin axis is also 
parallel to the body Y axis, and (3) the pendulous type nutation 
dampers which are connected to the primary part and are constrained 
to move in a plane that is perpendicular to the nominal spin axis and 
at a distance Z, above the transverse inertial plane. It is assumed 

\-tJJ " V 

the dampers are hinged or pivoted about a torsion wire support which 
offers a restoring (spring) torque in addition to the dissipative 

7 

torque . 

A. Equations of Motion for the Case of a Constant Speed Rotor 

The system to be studied is shown schematically in Fig. lb. The 
orthogonal coordinate axes X,Y (spin axis), Z are fixed in the parent 
body with the origin at its equilibrium position center of gravity. 

The main body has mass M, polar moment of inertia B, and transverse 
moments of inertia A and C. Four, small pendulous dampers each having 
mass, m, are attached to the main body, and are free to move in a plane 
y - Z, perpendicular to the polar axis. (The position of two of these 
are shown in Fig. lb. The other pair would be aligned, in equilibrium. 
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along the + X axes). The angles swept out by. the pendulous dampers are 

The analysis shall initially be undertaken with, all four end masses 

being able to move, n = 1,2,3 or 4. Later we shall simplify our 

analysis by assuming n = 1 or <j> = <j> = 4> =*0 for the actual case of 

2 3 4 - 

SAS-A. Labelling the weights with subscripts, 1-4, the coordinates of 
the end masses can be expressed: 

z x = r Q + cos z 2 = - (r c + r 2 cos <j> 2 ) , z 3 = - r 2 sin <j> 3 , 

z ^ sin <p 4 , x-^ = r^ sin $]_/ x 2 = - r l sxn $2' x 3 ~ r o + r l cos ^3 

x 4 = - Cr Q + r 1 ) cos <j> 4 , y= Z .... . (1) 

All external forces except gravity, are neglected and even with that we 
shall neglect the effects of the inhomogeneity of the earth's field. 

The motion of the center of gravity of the system and motion about the 
center of gravity can then be separated. 

If cj x , ojy, and gj z are the angular velocities of the main body about 
the X,Y, Z axes respectively, and r£ is the position vector of the ith 
weight, the kinetic energy of the main body and the damper can be 


expressed as : 


.1,3 


_ J- . 2 + x 

Ao)^ + — 




1 n 
i E 


, Bok" + 'r Cu)2 + — M v M/cm + 9 v 

2 y 2 z 2 2 1=1 1 . C2) 


m. v 


From the definition of the system center of mass: 


'cm/o 


n _ 

m ^i = 1 r i/o 

M + l m. 

i = 1 1 


(3) 
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Where point o is the center of the coordinate system. The velocity 
of the various components relative to the system center of mass may 
be expanded: 

v = v +v< (4) 

m. . ra. . o/cm v ' 

l/cm l/o 

v M/cm = v K /o + v o/cm * * 

The components appearing in Eqs. (4) and (5) can further be expressed 
as: 

V ’ ^ • ' • • ' - ’ ' • • ; (6> 

^o/M = 0 (7) 

- - A m Ir* -/Av 

o/cm cm/o cm/o M + Zm 

Upon substitution of Eqs. (6) - (8) into Eq. (2), the kinetic energy 
of the main body may be expressed as: 

2T M = Aw 2 + Boj 2 + Cuj 2 + m (E r. . r . ) - m 2 (E r. . Er.) . (9) 

X y 2 ^ 1 ^ ^ 

M 

where M = M + Em, 
i 1 

The kinetic energy of the rotor is given by: 

T Rotor ~ j oj 2 + I R y + s) 2 + I R Z to 2 J (10) 

Where I R , i = x,y,z are the rotor principal moments of inertia and 

i 


- 6 - 



the rotor is assumed to be spinning about the Y-axis with a relative 
angular velocity s, with respect to the main spacecraft. 

If it is remembered that, e.g. the X-component of r^ is 

• 

x. + o).,z. - to_y^ , and use is made of Eqs. (1) and (2), the kinetic 

energy takes the rather involved form given in Appendix A. The 

equations of motion for this case can thus be expressed in terms of the 

quasi-coordinates (co ,oj ,u) ) and the angle swept out (<J>.) by the 

x y z 

8 

pendulous dampers according to: 


d / 3T . 03 3T , u> 3T 

* 'it 1 ' z 5S- + ' , sr“ I * 

x y z 


ft* A? ‘ Jjj""-- 


d ( 3T j 3T . 


dt 3w 


3(i) 


♦ « |2- - K 

z 3d) < 


( 11 ) 


3 _ ( _ 5 £, -«■ 52 

dt 3w y 3o) 

z 


(i) 3T — 

+ x = L 

3(i) z 

x y 


and 


( 

dt 


3T 


) - 


3T 
3<j) . 

l 


3 F 


3 V 
3 <J) . 


= Q.' 
x 


1 n 


a^i 3 * a 

where L^, i = x,y,z are the applied torques and n, represents the total 

number of the dampers; the viscous forces on the pendulous end masses, 

which vary linearly with the angular velocity, 4 k , can be derived from the 

3 9 

Rayleigh dissipation function, F; in this case, ' 


n 

2 i=l 


F = I i (k.<p 2 ) (12) 

i i 
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The restoring torque on each damper provided by the torsion wire support 
is represented by = - 8V where the potential energy, V, is 

' 3?7 

proportional to, the square of the angular displacement from the 
equilibrium position, 

V = I ? K.<J> 2 (13) 

2 i-1 1 i 

Eqs. (11) may be expanded and simplified with small angle assumptions 

made relative to the magnitudes of the $j_/ and with further assumptions 

that ^x/^y, “z/^y and ^i/ w y are small in order to obtain first order 

damping contributions, for the case of n = 2 dampers and for the case 

3 

of a constant speed rotor: 

^ wm \ 

A(I) X + (C - B) u)y(i) z - io z ^'RyS - 2mr-j£ w y (^ - i}^) ~ L x • • • (14) 

Bo)^ + (A - C) w z m x + mr^ (r Q + r^) (<J>^ + $^) = L (15) 

Cu z + (B - A)a) x u) y - mr-j^^ - <{> 2 ) + mr^C^ - cf> 2 ) co 2 

+ lR y w x s = L z (16) 

mr 2 (1 - m v A + m 2 r? ~ • 

1 ~ ' < * ) , : — — ~ mr,-tio + mr, (r + r-,)u) 

1 - Iz lo ly 

" M 

+ mri (r o + } ^ u 2 _ ™ 2 rJ ^ 2 + mr^^y 

M M 

= - (17) 
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mr^ (1 - — )<K + m ^ r ? ,j, + rnr + mr. (r + r, )id. 

1 - 2 —5- 1 z 1 o 1 1 

M M 


+ • mr, (r + ”^1 \ <b co 

1 o - >2 v 

M * 


2 _ 


m 

M 


2d w 2 - mr,-£w u) 

r l 1 y 1 y > 


- k 4> - K$, 


(18) 


where the total inertia terms have main body, rotor and damper components: 


A = A + I!r + I 

x x 


B - B + t r + x d 


C - C + *R + I d 

z z 

which can be expanded to yield: 


_ 2 2 
A = A + I R + 4m£ + 2m (r + r-^) 

x 


B = B + I R + 4m(r Q + r 1 )' 


— 2 2 

C = C + I D + 4m£ + 2m(r rt + r ) 

K o jl . 

z 


and £ _ 

M 


k is the damping rate constant 

K is the restoring spring constant of the torsion wire support 




I 

for Variable Speed Rotor 

--- ■■ ... ■ — ^ __ 

The Eulerian approach is considered here so that it will be 
possible to distinguish between reaction torques, applied torques, 
damping torques and external torques. The angular momentum of the 
main body can be expressed as: 

L m = (A + I d ) w x i + (B + l d }m y j + CC + I d ]u z k ..... C19) 

x y z 

whereas for the rotor: 

L R = r R U R 1 + X R W R 3 + r R k (20) 

-xx y y z z 

'Using the familiar relationship, 

dL - = dL . + wXL = IN 

dt I space dt ' Body 

the main body equations can be developed as: 

(A + I d ) w x + (0 u z CC + I d 1 - uy^CB + I d ) = N d + R x . . (21) 

x y Z y X 

(B + I d ) « y + “ z “> x CA + I d ) - (C + I d ) = N d + R + \ (22) 

y x z y y 

i • 

(C + X d >■ " z .+ w x w y CB + r d ) - o, x o) y (A + I d ) = N d + R z . . (23) 

Z y X Z 

and the rotor equations, similarly: 

X R v ®R + X R “Ir “z^R = N + r 

x 2 z y y x x..., (24) 
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J R “R + “z^R Z R " U x“R J R = N R + \ 
y y xx zz Yy 


( 25 ) 


X R W R + «X“R I R - “y^R Ir = N R + Rr 


(26) 


z z 


y y 


X X 


where I R , are the principal moments of inertia of rotor, (i = x,y,z) 


I, , are the moments of inertia of dampers about X,Y,Z 
d i 

N, , are the damping torques on the main body 
a . 
l 

R^ , are the reactions (torque) of rotor on main body 


N R , are the (applied) torques acting on the rotor 


Rj^ , are the reactions (torque) of main body on rotor 


L m , is the motor torque 


In Eqs. (21) - (26) the effect of all external torques (e.g. aero- 
dynamic, solar pressure, gravity-gradient) have been neglected. 
Substituting first order expressions for previously derived, 
Eqs. (21) - (23) can be written as: 


• i 


(A + I ) u + ai u ( C + I, ) - a) a) (B + I ) = 2mr £(<p - <b )co 

d x x y z d z y z d y 1 1 2 y 


+ R 


(27) 
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( 28 ) 


(B + I d )» y + 03 z w x (A + I d ) - W X 0) Z (C + I d ) 


= - n ( r o + r i) (<t> 1 + <l> 2 ) + Ry + Ln, 


(C + x d >b x + Vy (B + J d 5 " VV (A + T d ) 

z y x 


= mr^£.(<|>^ - 4*2) ~ - <£>2)^ + 


(29) 


The assumptions are, that reaction torques are equal and opposite, 

, and the inertia (motor) torques about axes other than the 
i 

spin axis are zero, i.e. N R = N R =0, but ideally, N R = - L 

x z y y 

i.e. the inertia torque about the spin axis is equal and opposite to 
the motor torque. . * • :■■*■*.*■* 

Now combining, (24) with (27), (25) with (28) and (26) with (29) and 
noting that. 


= u t s , oj r =m , and 0 ^ = 


= 


we obtain. 


Aw x + (C - B)ayu z - w z I R s = 2mr ^(^ - 4> 2 )w y 

y 


Bm + (A - C)u) z to x + I R s = - mr 1 (r o + r.^) (<j> 1 + 4> 2 ) 

y ° 


. . . (30) 
* 

. . . (31) 


Cw z + (B - A)o) x a) y + m x I R s = mr 1 £(({) 1 - <£> 2 ) -mr-^^ 

y 

where A = A + 1^ + I R , and similarly for B and C. 

. x x 


4>_)u 2 .(32) 

2 y 
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In addition, the general torque equation for the symmetrical 
rotor, Ci.e. 1^ 1^ 1 is obtained from Eq. (251, . 


1 R t u y + s) “ Lr 
Y Y 


C33) 


where the following approximation has been made: 

L m * L R >>_ *y 

y y 

When s = 0, then Eqs. .(30) - (32) are identical to Eqs. (14) - (16) . 


C. Stability Analysis for the Special Case of an Undamped Symmetrical 
Satellite 

We will consider the system without damping, and a spacecraft 
symmetrical about the nominal spin axis, i.e. I, = 0, C = A, and 

l 

I R = I R . Under these assumptions Eqs. (30) - (32) take the follow- 
x z 
ing form: 

A'w + (A' - B')u> £0 - w I_ s = 0 (34) 

x y z z R 

y 

B'u + I R s = 0 . (35) 

y 

A'o) z + a) x ajy(B' - A') + u x I R s=0 (36) 

y 

The general torque equation, Eq. (33) is again: 

X R (“y + k) " L R 

y y 


It should be noted that Eqs. (34) - (36) are exact with no restrictions 
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on the magnitudes of to and a> . 

X z 

After multiplying Eq. (34) by cj x and Eq. (36) by to z and adding then 
a first integral results as: - 

w 2 + to 2 = K = Constant (37) 

x z 1 

Eq. (37) indicates that the amplitude of the vector sum of the trans- 
verse angular velocity components is bounded during spin-up. 

Eq. (35) can be integrated directly, to yield: 


_ I R s + B'u (0) + T R s (0) 

y y y 

(0 = * - : - 

y b' 


(38) 


Upon substitution of (1)^ from Eq. (35) into Eq. (33) , 


s - 


\ 


= c = constant . 


(39) 


* 1 - _U 

y B' 


Upon integration of (39), 


s = ct + s (0) . (40) 

■v 

which states that the relative angular velocity of the rotor, for the 
case of a symmetrical satellite increases uniformly during spin-up. 


Substituting the values of and s into Eqs. (34) and (36) we 
can get. 
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(»' - B') [»„(0)B' t t S(0)1 
*■“* + IT ». 

X''’ 

+ [ — ll2Jl T V - J V fct + S(0)J U =0.(41) 
B’ 2 

A « m + (B 1 - A*) [W v (0)B ( + Z R S(0)] 

A w z + 1 1 

B 1 X 

f (A* - B«) ( _I Ry) _ Z Ry] [ct +s(0)] j x = 0 • (42) 

B’ 

If p = (A* - B’) [6J y (0) B 1 + I R y S (0) ] 

— - 

and, Q = (A' - B 1 ) (- X R y ) x 

R ’ R 

b y 

then Eqs. (41) and (42) can be written: 

A 1 u) + u> [P + Q(ct + s(0)}] =0 (43) 

A Z 

A'<L - w rp + Q(ct + s (0) }] =0 (44) 

2* A 

From Eq. (44) and after differentiation with respect to time, 

A' [P + Q{ct + s{0)}]u) z - w z A*cQ 

w x = [P + ;Q(ct + s (0) } ] z (45) 

After substituting Eq. (45) into Eq. (43), there results 
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’>■ If *' 


cQ _ P + Q { ct + s CO ) } , ^ 

W z JP + Q {ct + s C01 }J z * A' J . “z “ 0 *-46) 

Similarly by expressing a) z = f (tu x ) and differentiating, 

•• • , i 2 

w _ cQ ' b) + . P + Q (ct + s (0) K oj =0 (47) 

IP + Q { ct + s (0) }] 1 A' X 

Eqs. (46) and (47) are second order differential equations with 

variable coefficients and have a regular singularity at time t = t Q » where 

t Q ^_-£F=“ j S(0) Q + P j . The solution of this type of differential 

9 Q c 

equation can be represented by a power series expanded around the 

10 

regular singular point.- At least one of the solutions tc Eq4 
(47) can be represented as: 

.-A V ‘ _ 

w x = l a k (t + a) r + k (48) 

k = 0 

where 'a 1 is the regular singular point at time, 

t = - f s -(P^ 2— + , Hi 

Qc 

oo r + k — 1 

Then, u> x = E a^ (r + k) (t + a) (49) 

k = 0 

•• oq 

and o) = Z a^ (r + k) (r + k - 1) (t + a) r + ~ 2 ... . (50) 

X K 

k = 0 

Now Eq. (47) can be written as: 

w x 0 2 2 2 

Wx - cmr + ct + a) u x = 0 (51) 

Substituting Eqs. (48) - (50) into Eq. (51), there results, 
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£ a Cr + kl (r + k - 1} Ct + al 
k = 0 k 


r + k - 2 


CO 

I 

k = 0 


a^ Cr + k). Ct + a) 


r + k - 2 


n2 r 2 m 

+ £-£_ £ a k Ct + a) 
A' 2 k = 0 


r + k + 2 


= 0 


(52) 


The indicial equation is obtained by equating to zero the coef f iciencts 

of the terms in the smallest power in Ct + aj and assuming a f 0. 

o 

In this case the smallest power of Ct + a) , obtained by setting 
k = 0 is (t + a) r 


Then, 

a Ir Cr - 1) - rj » 0 • . C53) 

o 

Since it has already been assumed that a f 0, then the roots of the 

o 

indicial equation are: 

r = 0 , and r = 2 

. r ” 1 

Equating the coefficient of Ct + a) to zero, we obtain: 

( r + l)r - a^r + 1) = 0 C54) 

From consideration of Eq. (54) it can be concluded that a^ = 0 for 
both r = 0 and r = 2 . 

By equating the coefficients of (t + a) r to zero it is seen that 
i s arbitrary for r = 0 and &2 is zero for r = 2. Similarly by 

• r -f 1 

equating the coefficients of (t + a) to zero we get a^ - 0 

for both the values r = 0 and r = 2 . 
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IT + 1c 

; Now by equating the coefficient of (t + a) to zero, the 

following relation is obtained: . ' -*• 

a k + 2 (r + k + 2) (r + k + 1) - a^ + 2 (r + k + 2) 


/ 


r? 2 

, Q c 
+ ~rnr a 


k - 2 


= 0 


(55) 


.From Eq. (55) , the general recurrence relation between the coefficients 
can be derived: 

a k - 2 


Q 2 c 2 


k + 2 


(r + k) (r + k + 2) 


(56) 


For a selected arbitrarily (a ^ 0) and for the root, r = 0, the 

o o . 

other coefficients may be expressed as follows: 


a., = 0; a 0 = 0; a = 0; a A = - O^c 


2^2 a. 


a r = 0; 


A ,z 2.4 


.2_2 a 


a A = - Q c 2 ; a 7 = 0 ; a R _ Q~*c 


i4„4 




4.6 


A 77 2. 4. 6. 8 


; a 9 = 0; 


4 4 

a io “ 2^ 


4.6.8.10 


6 6 


; a 1]L = 0; a., 0 = _ Q c 


12 


A ,b 2.4.6.8.10.12 


etc. 


, a m into Eq. (48) the solution 


After substituting a Q , a^, a 2 , . . . . 
corresponding to the root, r = 0, takes the following form: 


“x = % + a 2 (t + a) - Q 2 c 2 

1 ~T rz 

2_2 


a Q (t + a) 


- 2^ 


2.4 

a 2 


A 77 4.6 


(t + a) 
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( 57 ) 


Q 4 c 4 a o 

A' 4 2. 4. 6. 8 


(t + a) 


8 


Q 4 c 4 

A' 4 4.6.8.10 


2 (t + a) 10 


Q u c u Q o 

A ,b 2.4.6.8.10.12 


(t + a) 


12 


Similarly, for r = 2 and selecting 3- Q ^ 0, with the aid of Eq. 

(56), the coefficients can be related as: 


2 2 

a-, =0; a = 0; a. = 0; a = - 

1 2 3 4 


o 


A* ^ 4.6 

3c - 0; a fi =0; a_ = 0; a fi = — ^2 

5 6 7 8 4.6.8.10 


a 9 = ° ; . a 10 = 0; >11 = ° ; a l2 


Q 6 c 6 a o . 

A' b 4.6.8.10.12.14 ’ 


e‘tc . 


Again after substitution of the above values of a Q , a^, a 2 , 


a into Eq. (48), the solution corresponding to r = 2, becomes: 

m 


03 = 


a Q (t + a)' 


Q 2 c 2 
’ A 72 


4.6 


(t + a) 




+ b (t + a) 10 - K 5L ^ ( t + a 3 

A' 4 4.6.8.10 { A ,b 4.6.8.10.12.14 


a 6 c 6 


14 


+ . . . 


(58) 


The complete solution of Eq. (47) is, in general, an arbitrary 
linear combination of Eqs. (57) and (58): 


03 = C-,0) + c„0) 

x lx, 2 x 
1 2 


(59) 
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Eq. (59) can be expressed, using Eqs. (57) and (58), as: 


to = C, [a + a (t + a) 

v 1 O 


(t + a) 4 - 

TZ n a jiZ 


4 6 


(t + a) 


+ W 


A' 4 2. 4. 6. 8 


(t + a) 8 - Q 4 c 4 


A 71 * - 4.6.8.10 


(t + a) 


10 


_ q!£ 6 


a o 


A ,b 2.4.6.8.10.12 


(t + a) 12 -. . . .) 


+ C [a (t-+ a) - 
2 o 


2 -fz <t + a) 6 + 2^4 _!2 (t + a) 

A' 2 4.6 A' 4 4.6.8.10 


10 


-.6 6 a ^ 14 

2- . 9 2 (t + a) + . . .] 

A ,b 4.6.8.10.12.14 


(60) 


for - oo < t < co 


It should be noticed that the solution converges for all finite 


values of time since there is no other singular point of the differential 


Eq. (47). 

Eq. (60) can also be represented as: 


(dx = K Cos [ 8£ . (t + __ a ) 2 + j 

1 A* 2 2 

where 

K = and K, = Tan' 1 f (C^ + C' 2 ) 

1 CosK_ 2 1 2A J 

2 Qc 

C' = f 2 and c' = El 

1 a - 2 C. 


(61) 
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It can also be seen that Eq. (61) can also be represented as: 


w = K Cos / M'dt 
x 1 o 


( 62 ) 


where 


M ' = [ct + s (0) ] 


Q 


Now if we assume that the solution of Eq. (46) has the form: 


w 2 = £ b. (t + a) 
k = 0 


r + k 


then by applying the method of Frobenius as before we obtain the 
complete solution to Eq. (46) as: 

,2 


u z = K Sin [ ( b + 3. ) + k ] 

1 A' 2 ' 2 

where the constants can again be related by: 


(63) 


K, Cos K Q£_ = c . + c ’ , ibL 

2A' 3 4 Qc 


Tan K 0 = C-,b 
Z 3 o 


where 


C* = ^2 

3 b 

o 


and C' = Ei 
4 C„ 


Eq. (63) can further be reduced as before to the form: 


w = K Sin / M'dt 
2 1 o 


(64) 


From consideration of Eqs. (62) and (64) it is apparent that: 


w* + ui£ = Kf 
x z i 


which compares directly with the first integral, (37). 


( 65 ) 
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When the nutation damper is physically constrained from movement, 
the magnitude of vector sum of the transverse angular velocity components 
remain constant during wheel spin up under the influence of a constant 
motor torque. The stability criteria for this system is based on the 
magnitude of the transverse components of the angular velocity. If the 
transverse components of the angular velocity do not exceed the initial 
values, then it is assumed that no serious stability problems would be 
encountered. The implication of Eq. (65) is the boundedness of this 
motion during spin-up. 

D. Analysis for Asymmetrical Main Body 

Assuming the system is undamped and has a symmetrical rotor, i.e. 

I =0 and I = I , but with A f C, Eos. (30) - (32) can he 
cl. R R 

x x z 

expressed as follows: 

A'w + (C' - B')o) w - u I n s = 0 (66) 

x y z z R 

Y 

B'w y + (A' - C'Ju z u x + I R s = 0 (67) 

y 

C'lu + (B' - A')u) w +u T s = 0 (68) 

z x y x R 

y 

and the torque equation for the symmetrical rotor is the same as before: 


\ i&y + S) - L r 

y y 


After multiplying Eq. (66) by u)„ and Eq. (68) by. w and adding there 

** z 


results : 


! J, 1% ■ A ' <“x w x> + c ' K w z> 

CI)„U_(A' - c* ) = 

X Z CO 


(69) 


Substituting Eq. (69) into Eq. (67) and after simplifying we obtain the 
result as: 


A' (w x u> x ) + B(a)yO)y) + C * <0J z tu z ) + L R = 0 


(70) 


If the torque control law during spin-up is as: Body + L R =0, then 

y 

Eq. (70) becomes: 

A' (od x J> x ) + C'(to z dd z ) = 0 . . . (71) 

which results. in: 

A'w 2 + C'w 2 = c = Constant .• (72) 

X ^ 

and the motion is bounded as in the case of a symmetrical main body. 

In general, the integration of (70) yields the following first 
integral : 

1 A'to^ + 0^2 + 1^ c * to 2 + L od„(u)du = Constant : . . (73) 

2 x 2 y 2 z ° y y 

E. Effects of the Nutation Damper 

Under the assmumptions of a symmetrical spacecraft and only one 


damper free to move, the equations of motion, (30) - (32), are as 
follows: 

Au> x + (A - B)tdytd z -ui z I R s = 2mrj£$ u (74) 

y 

+ I R s = - mr^ (r^ + r 1 ) (75) 
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I ■ . 

... ! . 

! 

/ 

A <I> + (B - A)io.,u),. + w I s = mr,£$ - mr£<|>, w 2 

z 'xy x R 1 i 1 1 y 

i y . . 

and from Eq. (33) , the torque equation is written: 

I R <i y + S) - L 

y y 

In addition the damper equation, from (17) is: - - 

2 •• . “ • 

mr 1 (1 - E ) 4*2^ + mr^r^ + r u y " mr^^ 

M 

* ” r 1 (r o ’ ) *?’ 

M 

+ mrj^uwOy = - k<j> y - K<j>^ 

Eq. (75) can be integrated directly, to yield: 

«v . « y , 0 > - : r » - V (0) 

fi B B 

“ltro + r l} 

- $ (o) . 

1 1 

B 

Upon substitution of from Eq. (75) into Eq. (77), there results: 

, , \ ♦ - * ± z °-+ ill i 

I I R % 1 

R (1-_X , B (1 ~ ~~ ) 

B B 


( 76 ) 


(77) 


(78) 


(7S) 


(BO) 
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The first term on the right hand side is constant, say, c, under the 
assumption that L„ i s constant. - ^ """ 


Upon integration of Eq. (80) 


mr, (r + r,) * 

S = ct + ° 4>1 + s (0) . 

""5 (1 - , 4 ^ ) 

B 


(81) 


It is seen from Eq. (80) that compared to the symmetrical case, the 
relative spin rate of the rotor in this case does not increase uni- 
formly; also that if the damper mass is zero (effect of damping neg- 
lected) Eq. (40) would follow directly. 


F. Analysis for the Case of a Variable Torque Law 

In this section an analysis is made assuming a variable rotor 
torque instead of a constant torque as described in the previous cases. 
The rotor torque is assumed to be a linear function of time and reaches 
twice the value of the constant torque case during spin-up. It is 
clear that the average torque for the variable torque case is the same 
as that for the constant torque case previously treated. If we assume 
a symmetrical satellite without the effect of the damping, the equa- 
tions of motion can be obtained from Eqs. (34) - (36) as: 

A*m x + (A* - B») uy) z - s = 0 . . . . (82) 


+ 



( 83 ) 


( 84 ) 


A’w + (B' - A') w x o)y + co I R s = 0 . . . 

z x y 

and from Eq. (33), the general torque equation: 

I R (co + s) ~L R (85) 

y y y 

In addition, for the variable torque law: 

L r = ct (86) 

y 

After multiplying Eq. (81) by co and Eq. (83) by co and adding, 

X z 

there again results the first integral as: 


co 2 + co 2 = Constant 


(87) 


(87) indicates that the amplitude of the vector sum of the transverse 
angular velocity components is bounded during spin-up, as obtained before 
Upon substitution of (Oy from Eq. (83) into Eq. (85) and using Eq. (86) , 


s = c't 

where c' is a constant equal to 

Upon integration of Eq. (88) , 


( 88 ) 


J R (1 - X R ) 

y . y 

B' 


s = ^ 2 + s (0) . . (89) 

from which it can be concluded that the relative angular velocity of 
the rotor increases parabolically with time. 
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G . Effects of the Misalignment of the Principal Axes from 

the Geometrical Axes 

It is noted that the rotor will be displaced from the 
system center of mass and that the rotor spin axis will 
nominally be parallel to the desired main body spin axis. 

In the event that a perfect static mass balance is not 
achieved, the center of mass of the system will not be at 
the geometrical center of symmetry when the system is in 
equilibrium. For this situation there may be a misalign- 
ment of the principal main body axes from the geometrical 
(x,y,z) axes. In addition, if the rotor is not perfectly 
mass balanced there may also be a misalignment of the rotor 
principal axes from the geometrical axes. 

To consider these possibilities the first order non- 
linear equations of motion were developed for both the main 
body and rotor, as before, but now including all the cross 
products of inertia terms for both the main body and rotor 
as referenced to the x,y,z system. The main body and rotor 
equations may be combined as in Section B and expanded to 
include first-order small amplitude nutation damping 
effects. The following first order nonlinear equations 
result in the coordinates: o) x , Uy , u> z , s, and 


2 7 



/ 


[I XX + X R + Am2 + 2m(r o + r l-> 2 ]“x " [ I xy + \ + r ^ J * 


a R Jt0 Z " a R- * 


' r mr i vr o +r l ; *-i Jw x M y 

XZ X 


[I zz + 4«i* - I R - I + X R - - 2m(r 0+ r^ 2 ] 


“y“ Z 


+ [l yz + I R ]a, y + t X R “ z + X R “y Is - t X y z + X r ]“* 
yz yy yz yz 


*[l X y + I R ]s) z ii) x + m£r^(j) 2 + 2mr^£iJ^t 
xy 1 


+ 2mr, (r + r. )i, to 
1 o ± 1 z 


[I + !r + m£r 1 4» 1 ]i x ' ! + [I y % + I R + 4m (r Q +r ) 2 ] J) y 

^ vv j UV 


" Cl yz + X R + !r * + mr 1 ( r 0 +r 1 ) ^1 


U xy + * R + m5 ' r 1 4> 1 la) y U, z ‘ [lxx + J R 


I ZZ ~ I R la) z W x “ ^ I yz + I R ^x w y 
zz yz 


[l xz + I R ] U x " “z 2 > " ^R u x.“ X R ^z^ 3 
xz yz xy 


+' 2mr 0 r 1 (j) 1 <j) 1 (i)y 
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- ['x! + X R + mr l <r o +r l ) +i 1 “x - [I yz + J R '“i 


X z 


yz 


+ tl zz + I R + W 2 +■■ 2ni(r o +r 1 ) 2 ](i z - I R s 

zz yz 


- Mjltj - [I + Ir !“ x “ 2 + U x + Ir ]» x 
* yz xy 


- [I + Ip - I - I_ - ] oj ui 

yv K xx R x 

” yy xx 


" tl xz + *r + ^ 1 (r o+ r 1 )<j> 1 ] 


xz 


w y w z 


- [I xy * X R + mr l^l ] “y ' !l R “y + X R “x ls 
xy * xy yy 


- 2 mS,r 1 4 ) 1 (J> 1 (i)y 


( 92 ) 


Torque Equation 

- J S “x + X R "y - J R “z + J R » 

xy yy yz yy 

- L r + E R - (I R - *R >u z“x + *R “y“z 
y y xx zz xy 


2 2 

+ It-, w s + It, w - I„ (u - I wo) 
Rz Rz R x Rxy 

xy xz xz n yz 


-I 0 ) s 
R X s 

yz 


. ( 93 ) 
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Damper Equation 


m r i ( r 0 + r i ) Wy - mr 


1 A«*» Z + mr 1 (l - ® )<j, 


p _ 

= - mr 1 (r Q + mr l )4>i w y - mr^o^oiy - k^ - K<j» 1< . (9*0 


M 

Eqs. (90)-(9U) for the general case where all (or 
most) of the cross product terms appear can not be readily- 
solved analytically and must be evaluated by numerical 
integration. However, these equations will be considered 
for three special cases which are representative of numerical 
cases studied in Part III for the SAS-A and OSO spacecraft. 


Case 1 


The system will be considered for the case where __ 
J xz =I zx J xy = = 0. C ' = A', I = I and all 

J Z Z XX 

effects of rotor axes misalignment will be neglected. When 
it is also assumed that u /oj v << 1, u> /to <<1, y<<A’, and 

a j Zi y 

y<<B', the equations can be written, for the case of no 
damping, as: 

X ( ft » __U t ^ 

u y^Z " w Z -‘•R - ' '-Z -X~y‘ 


A' 0 ) v + ( A ’ -B 1 ) u v u z - U> 7m I_ s = y ( w„ + w to ) . . . (95) 


yy 


B 1 ajy + Ip s = 0 

yy 


( 96 ) 
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A'd> z + (B' -A' )u) v tD, r + U) Y 1^' s = 7 (w v - Uiy0) z ) . . .(9T) 


x y 


x R 


yy 


J B <“y + s) - L r 

yy y 


(98) 


[ n ] 

xyz 


It should he recalled that A' = A + 1^ and similarly for 

xx 

B’. Case 1 is dynamically similar to that of a spacecraft 
having no misalignment in its principal axes hut a small 
difference between each of its moments of inertia about 
the two transverse principal axes, (see Case D, p. 22). 

Consider the matrix whose elements are the principal 
and cross products of inertia, for this case, 

A' 0 y 

0 B' 0 

Y 0 A' [ 

xyz 

It can he verified that the eigenvalues of this matrix are 
B', A' + y respectively and that the y (or b^) axis still 
remains a principal axis hut that the new principal axes 
in the transverse plane are now misaligned from the x,z 
geometrical axes. 

An approximate solution to Eqs. (95)-(98) has been 
obtained using pertubation techniques and will be developed 
later in this section. 
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Case 2 

The following are the approximate equatipns of motion 

where I = I = I = I = a. I„_ = 0. C ' =A ' , It? = 1^ , 
xy yx yz zy - *' zz ^ ■ 

no misalignment of rotor principal axes, and no damping as 

before: 


A '“x 

+ 

(A* 

-B')o) to z - u z I R s = o[u 




yy 

B' fi) 

y 

+ 

J R 

s = a[a) x + to z + to ( oj z • — to 



yy 

A ' “z 

+ 

(B' 

-A' )to x a)y + to x I R s = a [to. 




yy 

Z R 

( 

+ 

>> 

3 

s) = L„ ......... 

yy 



y 


( 100 ) 


(l02 ) 


Case 3 


If no misalignment of the main body principal axes from 


the geometrical axes is assumed, but that I =3 and also 

R yz 

R and I = I = 0, the approximate un- 
xx xy . R xz 


C ' =A * , I = I 


zz 


damped equations become: 


X! 

• 3 

<! 

+ 

(A'-B* )to y to z - 

M 

tsi 

3 

1 

= 

3 [ a> 2 + U)yS - 

■] 

B' to 

y 

+ 

! r s = B[u z 

- w x 



yy 


• 3 

< 

+ 

( B ' -A ' ) t0 x (0y 

+ to 

X 


• 

. 

# 


yy 


(103) 

(10U.) 


yy 


x^z 


I R (o) y + s ) = L r + 6[(» z - w x ( 0 y - to x s] (106.) 

yy y 

In all three cases above I is the same as I as before 

R yy R y 
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It should be noted that Cases 2 and 3 reflect situations 
where the polar axis of symmetry for the main body and 
rotor, respectively, are no longer principal axes. These 
three cases correspond to cases considered in the next 
section using numerical integration techniques. 

An approximate analytic solution can -be developed for 
Case 1 using pertubation techniques. Eqs. (96) and ( 98 ) in 
the variables uiy and s correspond to Eqs. (33) and (35) 
developed in the stability analysis of an undamped symmetrical 
satellite. The solutions can be represented as before: 

s = ct + s ( 0 ) (10T) 

where c = L R /I R [l - I R .., /B']r 

y yy yy 

and 

Wy = 0)y(0) - C 1 t . ( 108 ) 


and 


c = I c/B' 

R yy 


After substituting the solutions of Eqs. (107) and (108) into 
Eqs. (95) and (97) for the case where s(0) = 0, there results 


A'(d x + (P + Qct ) to = y [ u» z + o) x (w (0) - Cjt)] 


(109) 
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A'd>„ - ( P+Qct ) 0 ) = y [ d) - a) (u (0). - c,t)] .. 

z, a A z y 


(no) 


where P and Q have "been previously defined after Eq. (42). 

For the case where y << A* and y << | B * —A 1 | a solution 

using perturbation techniques may be developed by substituting 

the zeroth order solutions for ta and , Eqs. (6l) and 

X z 

(63) , (previously obtained for the homogeneous form of Eqs. 

(109) and (110) into the right hand side of Eqs. (109) and 

(110) . The zeroth order solutions can be represented by 

= K c o s ( x + K 2 ) (Ill) 


0 ), 


u z = K sin (x + K ) 
o * 


( 112 ) 


where x = x(t) = 


Qc 

A’ 


( t + a ) ‘ 
2 


Eq. (109) may be differentiated term by term with respect to 

time and ui appearing on the left side eliminated by using 
z 

Eq. (llO) to yield: 

A * 0 ) + [ P+Qct ] {co [P+Qct ] + [to - oj (to (0) - c , t ) ] } 

X A’ A' o 0 ^ 


+ Qcw = y{o) + [u> (0) - Cit]w T - c n c o Y } . . . . (113) 

^ z y - 1 - x x -*• 





J 

After elimination of to by using Eq. (109.) and algebraic 

z 

simiplification of all terms appearing with the coefficient 
"v" . not 3 ner t.tint. • 

x *= |qc|( t+a ) /A ' , 0) x = - to z x, to z = to x x 

0 0 o o 

and a = F/Qc, there results: 



T<Ax o r \ ( °) + ci a (HU) 
A' L t+a J 


Following the analogous procedure beginning with Eq. (110) 
and eliminating ui x and to x terms on the left side, a second 
order differential equation in u> z may he obtained: 




CO 


z 


<0 _ 2 2 

+ (Qc) ( t+a ; “ z = 

t + a A' 


^ Z Q r COy ( 0 ) + c j a 

A ' t + a 


(115) 


Eqs. (llU) and (115) are nonhomogeneous differential 

equations with variable coefficients. It is clear that the 

solutions to the homogeneous parts of (llU) and ( 11 5 ) are 

the same as developed previously, namely Eqs. (6l) and ( 63 ). 

The particular solution can be obtained by using the method 

of the variation of parameters. It is convenient to write 

the complementary solutions as 

to Y = K' cosx + K' sinx . . . . (ll6) 

h i 2 

<o„ = - K' cos x + K' sin x (117) 

h 2 1 

where K' ,= K cos K , K' = - K sin K 
1 1 2 2 1 2 
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The particular solution for 0 ) is assumed to have the form: 

x ' 


ai = u n (t) K' cos x + u (t) K' sin x 
X p 1 1 2 2 


(118) 


subject to the 

u 1 K' cos 
1 


constraint 

x + u_ K ’ 
2 . 2 


After differentiating Eq. 

to time the two equations 

u, and u with the result 
1 2 


that 

sin x = 0 . . . ( 119 ) 

( 118 ) term by term with respect 
may be solved simultaneously for 


f 

u = K 3 sin 2x + K 3 K 2 sin 2 x . (120) 

1 ITx 2K ' x 

1 

f /. -v 

* = - ^-3^1 cos x _ ^3 sin 2x (l2l) 

2 2K ' x kx 

2 


The integration of these equations may be facilitated by 
noting that , 

1/2 ..- 1/2 

dt = [A'/2 |Qc|] x dx (122) 


and performing the integration of the right side with respect 
to x instead of t . 

The integration of these functions is accomplished by 
using relationships in Section 2.632 of Ref . 11 (see Appendix 
D) which are valid for x>0 and a = P/Qc>0. These integrals 
involve products of exponential functions with incomplete 
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gamma functions. A comprehensive discussion of the incomplete 

gamma function is given in the text by the Bateman manuscript 

. 12 . . . - 

project where it can be noted that the incomplete gamma 

function with a complex argument, r(b,ix), can be related 
to the Boehmer integrals C(x,b) and S(x, b) according to: 

r(b,ix) = e i7Tb/2 [ C ( x , b ) - iS ( x ,b ) ] (123) 


and 


r(b,-ix) = e i7rb / 2 [c ( x ,b ) + iS(x,b)] . 


• • ( 124 ) 


where 


C ( x , b ) = r t b_1 cost dt . . . . (125) 

x 


S(x,b) = /°° t^ sint dt 


• « • 


( 126 ) 


(for the real part of b<l) and the Boehmer integrals may 
be evaluated by the following series:-*-^ 

, oo , . m 2m+b 

C ( x ,b ) = F(b) cos (^iL ) - E A ( * • • • • 

K 2 m=0 (2m)! (2m+b) 


S(x,b) = f(b) sin ( bjr_ ) - E ) 


m 2m+ 1 + b 


(127) 


( 128 ) 


2 m=0 (2m + 1)! (2m + 1 + b) 

It can be verified that for the integrals appearing in Eqs. 

(120) and (121), b = - i. and no imaginary terms will 

2 . 

appear in the final answer. After some algebra and 

simplification, utilizing Ref. 11, Eq. (123) and Eq. (124 ) s 
it can be shown that: 
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u l (t) = (K^/K') '{l(l)cos K 2 + .1(2). sin Kg'- I 3 sin Kg} (129) 

i 


/ \ 
u 2 ^ > 


f fw /tm \ r t / i \ 

Viv^/.rv^ ; i x vo. ; 


bill JV o “ 


-r / /\ \ 

j- v ^ y 


cub iv ^ - 


-r / ^ \ 

i V J ) 


rr 1 
XVg J 


t ~\ o n ^ 
\ J_ D w /. 


where 

= K 3 K 1 /A'/2|Qc| ./2; K 3 = y [ w ( 0 ). + c ^ ] /A ' . . . . . (131) 


1(1) = i [ S ( 2x , - I)| x ] (132) 

'2 2 x = 0 

t 

1(2) = - 1_ [C(2x, - _L )| ] . . (133) 

/£” 2 ; l x =o 

"C 

-1/2 x 

1(3) = - i- x | (134) 

4 x =0 

t . 


The complete approximate solution for w x can then be 
represented by : 


= K cos(x + K p ) + K' u, (t) cos 

1 


x + K 1 


u. 


(t) 


sm x 


(135) 


Following the same procedure as explained above the 
complete approximate solution for qj, .may be developed as: 

£j 


0 ), 


= K 


s m 


(x + 


V 


K' u (t ) cos x + K' u, (t) sin x .(136) 
2.3 14 


where 

u 3 ( t ) = -(Kj^/K^M- 1(1) sin K 2 + 1(2) cos K 2 - l(3) cos K 2 > 

.. (13T) 
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u^(t) = (K^/K^Hld) cos K 2 + 1 ( 2 ) sin K 2 + l(3) sin K 2 > 

The constants and can be related to the initial conditions 

• • '' ’ 

on o> x , u) z , oj x > and w z * 

Attempts to obtain complete approximate solutions using 
perturbation techniques for the equations of Cases 2 and 3 
have not been successful. For these cases the form of the 
"co " equation and/or the torque equations are more complicated 

«7 

than those appearing in Case 1. Approximate solutions for 
rotor spin rate, s, can be obtained again in terms of the 
imcomplete gamma functions. When this solution is substi- 
tuted back into the left side of the w x and o> z equations 
&.***.* farther integration or differentiation of this complicated 
form would be implied. 


39 


III. NUMERICAL RESULTS 


i 

Since the SAS-A dual-spin spacecraft system was used as the basis 

for the mathematical, model in the previous analyses , the fcl lowing 

SAS-A spacecraft design parameters were used in the numerical calcu- 
13 


lations : 

Rotor spin rate 
Main body spin rate 
Satellite mass 
Polar moment of inertia 
of main body 

Transverse moment of inertia 


Illy 

M 

B 


of main body 
The motor torque 


lU 


A =. C 


J R 


2000 r.p.m. or 209 rad/sec 
0.5 rad/sec 
132.33 kg 
28.5I4 kg-m^ 


27.00 kg-m 
0.8 oz-in. 


A.l Calculations of Spin Rate of Rotor During Spin-Up for the 
Symmetrical Satellite 
From Eq. ( 39 ), 


s = 


l R 


= c 


k R 


V (1- -gf) 

8 x 2 . h5 x 9.8066 

1^ x 2.2 x 11.519 (1 - 11.55159 ) 

28.551 59 


0.^9171863 rad/sec* 


- to - 



! 


and the time for spin-up. 


t = 


; s - s (0) 


c 


_ ' 209 - 0 


.49171863 


= 425 sec or 7.08 min. 


The total time of spin up is 7.08 min. which is compared to an orbital 
period of approximately 90 minutes. 

During spin-up for the symmetrical satellite without damping, 

Eq. (38) can be used to calculate the change of main body angular 
velocity about the spin axis. 


a) - (0 (0) =. - R y [s - s(0)j 

3' 


Substituting the values of 1 R , B* , and s we obtain. 


“y ~ = 


- 011519 


28.551519 


x 2000 


= 0.8068922 r.p.m. ' 

In the actual SAS-A Satellite, its spin rate was observed to be 5 r.p.m. 
immediately after launch. The wheel was then uncaged and accelerated. 

g 

This resulted in a decrease in satellite spin rate to about 4 r.p.m. 
This change -in observed spin rate during spin-up compares with the 
0.8068922 r.p.m. in this calculation. 
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.2 Humeri cal Results for 0S0 - Spacecraft 

In connection with the computer studies discussed later numerical 
results are calculated using the following 030 spacecraft design 
parameters : 


Rotor spin rate 

s 

6 r.p.m. or 0.628 rad/sec. 

Satellite mass 

M 

6*1.25465 lb 

The motor torque 

L «y 

2 ft - lb. 

Polar moment of inertia of 
main body 

B 

34 slug - ft^ 

Transverse moments of inertia 

A 

131 slug - ft^ 


C 

136 slug - ft 2 

Polar moment of inertia of 

J R 

378 slug - ft 2 


rotor M 
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(Eq. 39) can be used for approximating the spin rate of the rotor 
for this case. - 

Prom Eq. (39) , 


s = 


Ry 


= c 


T (1 - ) 

Ryy B' 


378 (1 - jgf > 


= 0. 064113 rad/sec/ 


The time for spin-up. 




s- s(0) 
c 


0.628 

0. 064113 


=9.8 seconds 


- 43 - 



B. Results of Numerical Integration- 

In this section the results of numerical integration of the non- 
linear differential equations of motion for the most general case, 
i.e. the asymmetrical main part and the symmetrical rotor and also the 
effect of damping are presented. The purposes of the numerical investi 
gation are twofold: first, to verify some of the previous analytic 

results and, secondly, to compare the motion for different cases 
considered. The numerical integration was carried out using the IBM 
1130 and IBM 360/50 electronic computers. The RKGS and SIMQ sub- 
routines are used to integrate five nonlinear equations, i.e. Eqs. 

(25), (30) -(32), and (78). It should be noted that Eq. (25) is the 
more general form of the rotor torque equation. The subroutine RKGS 
uses the fourth order Runge-Kutta method for the solution of initial 
value problems. The purpose of the Runge-Kutta method is to obtain 
an approximate solution of a system of first order differential 
equations with given initial values. It is a fourth-order integration 
procedure which is stable and self-starting? that is, only the 
functional values at a single previous point are required to obtain 
the functional values ahead. For this reason it is easy to change 
the step size at any step in the calculations. The entire input of 
the procedure is: (a) lower and upper bound of the integration 

interval, initial increment of the independent variable, upper bound 
of the local truncation error; (b) initial values of the dependent 
variables and weights for the local truncation errors in each 
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component of the dependent variables; (c) the number of differential 
equations in the system; (d) as external subroutine subprograms, the 
computation of the right hand side of the system of differential 
equations; for the flexibility in output, an - output subroutine. The . 
SIMQ subroutine is used to solve the simultaneous system equations, 
for the accelerations: u) v ,(o„,u) , s and cfuin terms of the angular 
velocity and position coordinates. A complete Fortran listing of the 
computer program is given in Appendix B. 

In all numerical results to be presented here, the main body is 
assumed to spin with an initial component of 0.5 rad/sec. and one of 
the components of the transverse angular velocity, i.e. uj x ( 0) is 
chosen initially to be 0.000159 rad/sec. All other initial variables 
are chosen zero. 

In the first case considered, the spin-up for the satellite with 
symmetrical main part without damping is shown. Fig. 2a shows the 
linearity of the rotor spin rate with respect to time. It also shows 
the rotor reaches its nominal spin rate of 209 rad/sec after a time 
interval of 425 seconds as previously calculated. Fig. 2b illustrates 
the time history of the transverse components of the main body angular 
velocity. It is seen that these components have a constant amplitude 
of 0.000159 rad/sec., the initial value of ijj x . Therefore the first 
integral expression, Eq. (37) has been verified. It should be noted that 
the time response of these components is that of a compressed 



sine wave with increasing frequency. This is explained from Eqs. (61) 

and (63), which show that the frequency increases directly with the 

square of the time. Since the components of the transverse angular 

velocity never exceed the initial value, from the concept of stability 

previously explained, no serious stability problem would be encountered 

here. It should be noted that at the completion of the wheel spin-up 

manuever, the transverse angular velocity magnitude of 0.000159 rad/sec 

could be removed by activating the nutation damper as described 
3 

previously. It can be seen from Fig. 2c that the main body spin rate 
decreases linearly with respect to time, which verifies Eq. (38). 

In the actual orbital configuration of the spacecraft a small 

asymmetry in the main body exists. Fig. 3 shows the motion of the 

* ■ ■ -.< * + . '? • 

transverse components of the angular velocity for the case of a small 
mass asymmetry in the main body without the effect of the damping. 

The moment of inertia parameters coorespond to SAS-A early design 
parameters. A small increase in amplitudes of both the transverse 
components of angular velocity is noted from the figure. No signifi- 
cant nonlinearity of the rotor spin rate nor the main body spin rate 
is observed from the data obtained by the computer simulation; there- 
fore these graphs are not shown. 

The total computer time for running the symmetrical case and for 
the case of small asymmetry varies from 20 to 25 minutes with the 
IBM 1130 computer. When large mass asymmetry in the main part was 
included a significant increase of IBM 1130 computer time was noted 



due to the excessive iterations required in the RKGS subroutines to 
maintain the same accuracy as in the previous two cases. So for 
better performance the numerical integration for this case was per- 
formed using the IBM 360/50 computer instead of the IBM 1130. 

Figs. 4 illustrate the effect of large mass asymmetry in the 
satellite main part without the effect of the damping. From Fig. 4a 
which shows the time history of the transverse components of the main 
body angular velocity, during the first 20 seconds, it is seen that one 
component of the transverse angular velocity has an amplitude almost 
twice the initial value. Since the stability criteria is based on the 
boundedness of the transverse components of the angular velocity, a 
problem of stability could be encountered in this case, especially 
in the presence of external torques which are continuously acting on 
the main spacecraft. A significant phase change in the time response 
of the transverse components of angular velocity when compared to the 
symmetrical case is also noted for the case of large asymmetry. Along 
with the phase change a significant increase in the frequency of the 
transverse components of the main body angular velocity is also noted. 
For the above two reasons the . response is shown only in the intervals 
from 0 to 20 seconds, and from 400 to 420 seconds (Fig. 4b), compared 
to the total interval from 0 to 450 seconds for the case of symmetrical 
satellite. For the rest of the time interval (not shown) the time 
response does not show any significant change in amplitude variation 


. from that shown. 



The effect of the nutation damper for the case of spin-up with 
a symmetrical main body is shown in Figs. 5. Fig. 5a is the repre- ^ 
sentation of the time history of the transverse components of the 
main body angular velocity. A small increase in the initial value 
of the amplitude is noted for both the components. No significant 
reduction in these initial amplitudes is noted due to the action of 
the nutation damper. Fig. 5b is the response of the nutation damper 
during spin-up. The figure shows a bias around the value 0.006 rad. 
The reason for the bias can be explained by the fact that during the 
derivation of the equations of motion the lateral center of mass 
shift due to the damper motion was not included. V7hen one damper 
is free to move, this shift could be more noticable than when^a pair - 
of dampers, diametrically opposed in equilibrium, is used. In the 
actual SAS-A post-launch performance a small damper bias angle was 
actually observed. 6 » ^ One of the causes for this phenomenon was 
the actual lateral shift in the spacecraft center of mass due to 
small errors in the final mass balancing prior to launch, It is 
also seen from Fig. 5b. that the nutation damper reaches a value of 
1.0168° at 425 seconds. 

Figs. 6 illustrate the damping effect considered for the case 
of small mass asymmetry in the main body. A small decay in amplitude 
of one of the transverse angular velocity components (w z ) is observed 
from Fig. 6a. If the initial rate of decay were extrapolated linearly 
the time constant for the decay would have been 31.79 min. (where the 
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time constant is the time required to reach 1/e of the initial 
amplitude). When compared with Fig. 5a, a small phase shift is 
observed. Fig. 6b, which is the damper response for this case, shows 
a noticable growth of damper angle amplitude during the time of spin- 
up. It can be concluded that the activation of the nutation damper 
during spin-up for this case results in an improvement in the system 
stability. 

When the effect of the nutation damper is included for the 
spacecraft with large asymmetry, it is seen from the Figs. 7a and 7b, 
that a small average decay of the amplitude of one of the transverse 
components (co x component) of main body angular velocity results. The 
amplitude of the iu 2 component never exceeds the amplitude during the 
first cycle of oj z motion. If the initial rate of decay were extra- 
polated linearly the time constant for the decay would have been 

41.06 min. This result could be compared with the time constant of 

17 

SAS-A during nominal performance, which is 22.3 mins. Comparing 
the two cases, i.e. cases of large asymmetry with damping and without 
damping, it is observed from the Figs. 4a and 4b, that the maximum 
value of is slightly higher than the value of oj„ in Figs. 7a and 
7b. So it can be hypothesized that, for the case of large asymmetry 
without damping the energy is being transferred into the transverse 
motion from another mode. It is observed from the Fig. 7c, that the 
damper reaches a maximum value of 0.713° at 15.6 secs. The results 
presented in the Figs. 4 and 7 were obtained using the IBM 360 
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computer, requiring about 130 minutes of running time for each spin-up 
case. 

All the numerical calculations in the above cases are based on 
the assumption of a constant rotor torque of 0.8 oz-in. Two cases 
are considered where the torque is assumed to be a linear function of 
time and reaches a terminal value of twice the average value of the .. 
constant torque, (i.e. 1.6 oz-in) during spin-up. Figs. 8 and 9 show 
the effect of this variable torque law during spin-up. From Fig. 8a, 
which is the time history of the transverse components of the main 
body angular velocity for the case of a symmetrical satellite without 
damping, it could be concluded that the vector sum of the transverse 
components of the angular velocity maintains the same amplititude through 
out the motion as previously shown analytically, i.e. Eq. (87). 

Fig. 8b shows that the rotor spin rate is a parabolic function of 
time and reaches the nominal spin rate of 209 rad/sec after a time 
interval 425 seconds as before. The main body spin rate also exhibits 
a parabolic variation as seen from the Fig. 8c. Under the consideration 
of the variable torque law applied for the case of a symmetrical satel- 
lite including the effect of the damping, a small decay in amplitude 
of one of the transverse components (w^) of angular velocity is noted 
from Fig. 9a. The time constant for the decay is calculated to be 
27.82 minutes. In addition a phase shift is noted when compared with 
Fig. 8a. In the case represented by Figs. 9b and 9c a small departure 
from the parabolic variation shown in Figs. 8b and 8c is noted from 
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consideration of the' printout , Because of the small magnitudes of this 
departure, differences "between Figs. 8h and 9t> and also between Figs. 

8c and 9c are not apparent within the plotting accuracy. For this case, 
the nutation damper reaches a maximum value of 1.356° at 425 secs., 
which is slightly higher than the maximum angle attained by the nutation 
damper for the case of the symmetrical satellite with damping under the 
influence of a constant motor torque (Fig. 5b). 

B . 1 . Numerical Integration Results with SAS-A Spacecraft and Principal 

Axes Misalignment 

In all numerical results presented above, the effect of the mis- 
alignment of the principal axes from the geometrical axes of symmetry 
is neglected. A few cases are considered here including the above 
mentioned effects. All the numerical integration of these cases were 
performed using the IBM 360/50 electronic computer. 

In all of the numerical results to be presented for the SAS-A 
Satellite with principal axes misalignment the main body is assumed to 
spin with an initial component of 0,5 rad/sec., and results are obtained 
varying one of the components of transverse angular velocity i.e. 
oj^. CO); also all other initial variables are chosen zero as before. All 
the moments of inertia (including the cross-products) parameters are 

13 

based on the Small Astronomy Satellite (SAS-A) orbital configuration. 
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In the first case considered, w x Col. is chosen 0.000159 rad/sec. 

i .... 

Fig. 10a, shows the time history of the transverse components of the 
main body angular velocity. It is seen from the figure that during the 
first 100 seconds the amplitude of both the transverse components of main 
body angular velocity increases to approximately 7 times the initial value. 
After the first cycle of motion, both the components of transverse angular 
velocity show a decay in their amplitude. Since the stability criteria is 
based on the boundedness of the transverse components of angular velocity, 
a problem of stability could be encountered here, especially in the 
presence of external torques which are initially and continuously acting 
on the main spacecraft. From Fig. 10b and 10c which show the time 
response of the main body and rotor spin rate s respectively , it is seen 
that no significant nonlinearity of both of these motions is observed within 
the plotting accuracy. 

The effect of the nutation damper for the case of Fig. 10 is shown 
in Fig. 11. From Fig. 11a, it is observed that during the first cycle 
of the motion, the amplitudes of both transverse components of main body 
angular velocity increase much greater than the initial value. After the 
first cycle of motion of the transverse components, significant reduction 
in these amplitudes is observed due to the action of the nutation damper. 

The time constants for the decay are 1.1+5 and 5-^07 minutes for w x 
and oi respectively. A small phase change in the time response of the 
transverse components of angular velocity when compared to the undamped 
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case is also noted. Fig. lib is the response of the nutation damper 
during spin-up. The figure shows a bias around the value 0.02 rad. 
which is greater than the bias observed when the effect of the inertia 
cross products, are neglected (Fig. 5b). It is also seen that the 
nutation damper reaches a maximum value of 3.391° at 42 6 seconds. 

Figs. 12 illustrate the motion of the satellite with the initial 
transverse component increased by a factor of ten. Other input conditions 
are kept as before. From this figure a small percentage increase in the 
amplitudes of both the transverse components of angular velocity is 
noted. After the first cycle of motion a decay in amplitude of both 
the components is also observed, but the rate of decay is less than 
that for the case of Fig. 10a. , The effect of damping in this case is 
illustrated in Figs. 13. An ’.increase in amplitude of both the trans- 
verse components of main body angular velocity is observed from Fig. 13a 
similar to that shown in Fig. 12. A decay in amplitude of both the 
components of the transverse angular velocity vector is noted with a pro- 
jected time constant of 12.27 and 50.97 minutes for io x and Wy. 
respectively. A phase shift in the time response of the transverse 
components of angular velocity when compared to the undamped case (Fig. 

12) is also observed. It is seen from Fig. 13b that the nutation damper 
shows a steady growth in amplitude after a few initial cycles of 
motion and reaches a maximum value of 6.87° at 425 seconds which is 
about twice the value obtained compared to the case where io^O) is 
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chosen to he 1/10 the value for this case. 

The initial value of one of the transverse components of main 
body angular velocity (<o ) was once again increased, to 100. times the 
value in Fig. 10a and numerical integration of the equations of motion 
was again performed. From the time responses (Fig. 1*0 no significant 
increase in amplitude nor of the decay in amplitude in either of the 
transverse components is observed. It is seen from all the time 
response curves of the transverse components of main body angular velocity, 
that the frequency of the motion increases as co x (o) increases. 

When the effect of the nutation damper is considered in this 
case, a small decay in amplitude of one of the components of the trans- 
verse angular velocity (oj x ) vector is noted from Fig. 15a with a pro- 

.'r ’ 'v 

jected time constant of ^8 minutes. From Fig. 15b, it is noted that 
the maximum angular displacement of the nutation damper is Ul° which 
is a physical impossibility since in reality the damper would hit the 
mechanical stops at +_ 20° amplitude. As expected, the larger value 
of initial transverse velocity (i.e. system nutation angle) results in 
a higher amplitude nutation damper motion. Because of the larger 
amplitude damper motion, the damper bias angle previously referred to- 
is not apparent from the scale used in Fig. 15b. 

The average IBM 360/50 running time for the cases shown in 
Figs. 10-15 was about 15 minutes. 

B.2., Numerical Results Using 0S0 Spacecraft Parameters 


All the numerical calculations in the above cases are based on 


SAS-A design parameters. The OSO series of dual-spin spacecraft is of 
current interest to many investigators. The numerical integration of 
the general equations of motion is carried out with the OSO spacecraft 
design parameters vising thelBM 360/50 electronic computer. 

In all numerical results to be presented here for the OSO 
spacecraft, the main body is assumed to spin with an initial component 
of 0.5761 rad/sec as calculated previously, and results are obtained 
varying one of the components of the transverse angular velocity, i.e. 
u x ( 0). All other initial variables are chosen zero. 

In the first case considered, to (0) is chosen 0.003170 rad/sec, 

A 

which corresponds to an initial system nutation angle ( the angle 
between the total angular momentum vector and the nominal spin axis ) 
of 0.1 degree. From the time response curve (Fig. l6a) of the trans- 
verse components of main body angular velocity, it is noted that both 
the transverse components complete about one-half cycle of their 
motion during the total time of spin-up. A small increase in amplitude 
of both the components during the first cycle of the motion is also 
observed. Fig. l6b illustrates the motion of the main body spin rate. 
It is seen from the figure that the main body is essentially de-spun 
inertially after the total time of spin-up of 9.6 seconds. It is 
observed from Fig. l6c that the rotor reaches its nominal spin rate of 
0.628 rad/sec after a time interval of 9.8 seconds in an essentially 
linear manner. 
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Fig. 17 illustrates the effect of a greater initial nutation angle. 
In this case the value of w„(o) is increased to 0.03172 rad/sec which 

A 

corresponds to a nutation angle of 1°. It is observed that Fig. 17 is 
identically in phase with Fig. l6a except for the scale factor on the 
ordinate. 

When the initial system nutation angle is further increased by a 
factor of ten (Fig. 18) , the phasing of the response of u x and u> z is 
still identical to "that shown in Figs. l6a and 17. The frequency 
response of the 0S0 main body components appears to be less sensitive 
to changes in initial nutation angle than that for the SAS-A system, 
perhaps because for the 0S0 system the majority of the momentum as 
well as inertia contribution is associated with the rotor. 

In the previous 0S0 cases considered (Figs. 16-18) the geometrical 
axes of the main part and rotor were assumed to be the principal axes. 
The effect of the misalignment of the principal axes of both the rotor 
and the main part from the geometrical axes of symmetry is illustrated 
in Fig. 19. In this case an initial nutation angle is chosen to be 
the same as in the case of Fig. 17. No significant change in any of 
the motion is observed as compared to Fig. 17. 

Some results are now obtained with an increased initial value of 
main body spin rate and no misalignment of principal axes. In all 
cases mentioned below the main body is assumed to spin with an initial 
component of 6.985^3 rad/sec and results are obtained varying one of 
the components of the transverse angular velocity, i.e. w x (0). 
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All other initial conditions are chosen zero as "before.- - 

In the first case considered, w x (0) is chosen 0.03172 rad/sec., 
which corresponds to an initial system nutation angle of 0.21°. From 
Fig. 20a a small increase in amplitude of one of the transverse 
components of angular velocity (m z ) is noted. After the first cycle 
of the motion, a decay in amplitude of w z is also observed from the 
figure. The rotor spin rate curve (Fig. 20c) shows that the rotor 
reaches its nominal spin rate of 0.628 rad/sec after time interval 
of 9-8 seconds as previously calculated. No severe nonlinearity of 
main body spin nor of the rotor spin rate is observed from Figs. 20b 
and 20c. The above mentioned results could be compared qualitatively 
with the cases in Fig. 3, where initial amplification is observed 
for both the transverse components of main body angular velocity. 

Fig. 21 illustrates the effect of a greater initial nutation 
angle. In this case the value of oj x ( 0) is increased by a factor of ten 
From Fig. 21 a small decay in amplitude of one of the components (w z ) 
is observed after the initial amplification of both the components. 

It is also observed that the rate of decay is less than for the case 
in Fig. 20a. Figs. 21 and 20a are identically in-phase except for 
the scale factor in the ordinate. 

The nonlinearity of the main body spin rate and the rotor 
relative spin rate (Figs. 22b and 22c) is observed when the initial 
amplitude of co x is again increased by a factor of ten. This 
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nonlinearity (present to some extent in all 0S0 cases but observable 
here due to the large initial nutation angle) is attributed to the 
fact that for the 0S0 spacecraft the rotor polar (principal) axis is 
no longer an axis of rotor symmetry. No significant decay in 
amplitude is observed from the time response curve of the transverse 
components of main body angular velocity (Fig. 22a), though it 
differs in phase and amplitude with Figs. 20 and 21. 

The effect of the misalignment of the principal axes of the main 
body as well as the rotor from the geometrical axes of symmetry is 
illustrated in Fig. 23. The initial nutation angle is the same as in 
Fig. 20, i.e. 0.21 degree. The amplitude of both the transverse 
components of main body angular velocity has increased by a factor 
of almost 2. After the initial amplification a significant decay in 
the amplitudes of both the transverse components is also noted. The 
linearity of the main body and rotor spin rate responses (figures not 
shown) is not significantly affected by the inclusion of inertia 
cross products terms in the general equations of motion. It appears 
that the nonlinearity in oj and s responses is more sensitive to change 

v 

in the initial nutation angle than to the presence of small cross pro- 
ducts of inertia terms. 
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B , 3 Comparison between the Numerical Integration and' Approximate 

Solution for the Special Case Considered _ 

In this section a comparison is made between the results of the 
numerical integration and the approximate solution for the case which 
is dynamically similar to that of a spacecraft having no misalignment in 
its principal axes, but a small difference between each of its moments 
of inertia about the two transverse principal axes (Eqs. 95-98). To 
compare the approximate solutions (Eqs. 135 and 136) with the exact 
solutions (numerical integration), the numerical evaluation of the 
approximate solution was carried out using the IBM 360/50 electronic 
computer. The numerical evaluation of the approximate solution for 
,the. SAS-A spacecraft during spin- up was performed for a total time 
interval of 450 seconds with a time step of one second. Two cases 
are considered here varying the only inertia cross product considered 

i . e . I « 

xz 

In the first case considered the value of I v _ is chosen 

xz 

2 

0.1048 kg-m . Fig. 24 illustrates the comparison between the approxi- 
mate solution and the exact solution for this case. The solid line, 
which is the result of the numerical integration shows a small 
initial amplification in amplitude of one of the transverse components’ 
(w x ) of main body angular velocity. No initial amplification can be 
observed from the results of the approximate solution (shown by 
dotted lines), but a small change in phase compared to the results 
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of the numerical integration is observed from Fig. 2k. 

The effect of the increased value of I (0.5) .is illustrated 

xz 

in Fig. 25. The results of the numerical integration show a significant 
initial amplification of both the transverse components of main 
body angular velocity. No decay in amplitude of any of the velocity 
components is observed within the specified time interval. In this 
case the phase difference between the results of the approximate 
solution and the results of the numerical integration is more noticable 
as compared to the case considered in Fig. 2k. A small increase in 
initial amplification of both the velocity components is also observed 
from the computer print-out of the approximate solution which is not 
noticable within the plotting accuracy. It could be concluded at this 
point that the larger is the difference between the moments of inertia 
about the two transverse principal axes , the greater is the change in 
phase between the approximate and the exact solutions. Of course the 
approximate pertubation solution can no longer be expected to provide 
reasonable convergence when c(|l xz | = . 

The case considered in Fig. 2k could be compared with the case in 
Fig. 10a, where an initial amplification in amplitude of a factor of 
seven is observed for both the transverse components. By observing 
the cases considered in Figs. 2k and 25, it could be concluded that 
when the nominal spin axis (polar axis) is no longer a principal axis 
(I^y i- 0, Iy z f 0) the large amount of initial amplification of the 
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transverse components of main body angular velocity is the result. The 
amplification is more sensitive to a small misalignment between the L 
principal axis of spin and. the nominal spin (geometrical -yy) axis, 
than to small differences between the two transverse moments of inertia. 
The amplification can further be explained as, when the nominal spin 
axis is no longer the principal axis , the motor torque has a component 
perpendicular to the nominal spin vector. This perpendicular component 
could cause an excitation in the motion about the transverse axis due 
to the altering nature of the torque component, depending on the 
relative phasing. 

The total time of execution of the approximate solution for 
both the cases is less than the running time of the numerical integration. 
For example, the total time of execution of the approximate solution for 
the cases when I Xz = 0.1048 and I xz =0.5 are 12.33 and 20.19 minutes 
respectively compared to 32.83 and 26.31 minutes respectively for the 
numerical integration. By performing the ratio test of the series 
(Eqs. 127 and 128) associated with the approximate solution, for the 
numerical case considered, it is observed that a minimum of l6 terms 
is required for the series to converge. By evaluating the approximate 
solution considering the number of terms in the series to be 16 ,20_, 
and 27 respectively, no changes within five' decimal places are observed 
from the computer print-out. It is interesting to mention further that, 
the computer is unable to evaluate the factorial in the denominator 
of the series when the number of terms in the series exceeds 27, 
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because of its sensitivity towards generating a large real constant 
(maximum magnitude ~ 10*^). Stirling's approximation is the best 
recommended method for evaluating a large factorial, and could be used 
in cases where series convergence is improved by increasing the 
number of terms. 


IV. 


CONCLUSIONS 


As a result of the present analysis and numerical results the 

following conclusions can be made: 

(1) For the case of a symmetrical satellite with no damping, the 
magnitude of the vector sum of the transverse angular velocity 
components remains constant during wheel spin-up under the 
influence of a constant motor torque. 

(2) With a small mass asymmetry in the main body, and without damp- 
ing, behavior is similar to the symmetrical case, but a small 
increase in amplitude of one of the transverse components of 
angular velocity is noted in addition to a phase change. 

(3) For the case of large asymmetry in the main body, one component 
of the transverse angular velocity has an amplitude approximately 
twice that of the initial value. Stability problems could 
result for this case in the presence of all the external torques 
which are continuously acting on the main spacecraft. 

(4) The effect of the nutation damper during spin-up is significant 
only for the case of an asymmetry in the main spacecraft, where 
a small decay in the amplitude of the transverse angular 
velocity vector is noted. There appears to be little advantage 
(or disadvantage) in activating the nutation damper for the case 
of no asymmetry. 
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For the case where motor torque is proportional to time and the 
spacecraft is symmetrical, the activation of the nutation damper 
during spin-up results in a small decay in the amplitude of the 
transverse angular velocity vector. 

When the effect of the misalignment of the main spacecraft principal 
axis from the geometrical (polar) axis of symmetry is considered, 
a problem of stability could arise due to the large initial amplifi- 
cation the system nutation angle. 

For the case of a dual-spin spacecraft with a large asymmetrical 
rotor, a nonlinearity of the main body and rotor spin time responses 
can result, depending on the initial nutation angle. This could 
cause an error in reaching the nominal terminal conditions. 
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v SPIN AXIS 



Figure la. Elements of SAS-A Attitude Control System 


(rotor displaced from spin axis) 
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Figure 2a. Rotor Spin Rate During Spin-Up 
Symmetrical Satellite and No Dampi 
. A = C = 27.00 kg-m2; B = 28.54 kg- 



•oas/avu 

‘.5 ‘miOH 30 Ail 3013 A RVinOMV 3AI1V13H 

- 69 - 


TIME, t,. SECONDS 



SQNO 





TIME, t, SECONDS 



Figure 4a. Time History of the Transverse Components of Main Body Angular Velocity 
Satellite with Large Asymmetry in the Main Body and No Damping 
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Figure 5b. Nutation Damper Response During Spin-Up 
Symmetrical Satellite with Damping 
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Components of Main Body Angular 
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Figure 7c. Nutation Damper Response During Spin-Up 
Satellite with Large Asymmetry in the Main Body and Damping 

A = 15.00 Kg -m2 
B = 28.00 Kg-m 2 
C = 2 . 00 Kg-m 2 
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Figure 8c. Main Body Spin Rate During Spin-Up 
Symmetrical Satellite and No Damping (Variable Torque Lav;) 

A = C = 27.00 Kg-m 2 
B = 28.54 Kg -m2 
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Spin Rate During Spin-Up (SAS-A No Damping) 

27.015; Iyy = 28. 5 1 *; I Z z = 26 . 905 ; Ixy = 2,123 x 10” 
0.10U8; I yz = 1.193 x 10-3; Ux (o) = 0.000159 rad/sec 












Fig. lib. Nutation Damper Response During Spin-Up (SAS-A with Damping 
1 xx = 27.015; Iyy = 28. 5^; Izz';= 26 . 905 ; Ixy= 2.123 x 1C“3 

Ixz = 0.10^8 Iy z = 1.193 x 10 -3 

*. m = 0.2158 kg; lc - 7.05 x 10- 3 ^ K = 6.10 x 10~ 5 0 ^( 0 ) = 0.000159 rad/sec. 
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Fig. 13b. Nutation Damper Response During Spin^Up ((SAS-A with Damping) 

X xx = 27 - 015 '> I yy - 2 8-54; I zz = 26. 905 ; 1^ = 2.123 x 10-3; 

I xz = 0.1048; I yz = 1.193 x 10 -3 

m = 0.2158 kg, k = 7.05 x 10“5, K = 6.10 x 10 -5 , w (o) = 0.00159 rad/sec 
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Time History of the Transverse Components of Main Body Angular Velocity 

(030 Spacecraft) . * ‘ . 

v 030 Spacecraft with Inertih Cross-Products 
I X y - 0.1; I xz - 0.0; Iy Z ” 

■Ip = Ip = 0.0; Ip = 0.0593 
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Fig. 20b. Main Body Spin Rate During Spin-Up (0S0 Spacecraft) 

■kji — 131.00; • - 31.00; I Z2 — 136. 00 ; Ijj xx 199*00 

% = 3T8.00 ; Ir zz = 208.00; ai x (0) = 0.03172 rad/sec 
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Fig. 20c, Rotor Spin Rato During Spin-iiUp (0S0 Spacecraft) 

; * I xx = 131.00; Iyy = 34.00; l| z = 136.00;. I Rxx = 199.0 

in = 378.00; I» = 208.001; uy(0) = 0.03172 rad/sec 






Fig. 21. Time History of Transverse Components of Main Body Angular Velocity 
. , (0S0 Spacecraft) 

I xx ~ 131.00; Iyy = 3U.00; I zz = 136.00 

I R XX = 1-99 • 00 ; Ir^ = 378.00; Ir z * 208.00 
to x (0) = 0.3172 rad/sec. 
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Fig. 22a. Time History of the Transverse Components of Main Body Angular 
! * Velocity (080 Spacecraft) . 

I xx = 131.00; Iyy = 34.00; l| z .= 136.00; I R =199.00 
iRyy = 378.-00; Ir zz = 208.00, mx(0) = 3.172 rad/sec. 
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Fig. 221 ). Main Body Spin Rate During Spirit Up (0S0 Spacecraft) 

^xx - 131.00; lyy - 3^.00; I zz = 136.00; 1 r xx = 199*00; = 378. 

I H z z = 208.00, u x (0) = 3.172 rad/sec. 
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Fig. 22c, Rotor Spin Rate During Spin- Up (OSQ Spacecraft) ■ 

Ixx = 131.00; Iyy = 3U.00; I zz = 136.00 ; . IR XX = 199-00 
lR yy = 378.00 Ir zz = 208.00, w x (0) = 3.172 rad/sec. 
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Fig 1 , 23*. Time History of the Trans vei^e Components of Main Body Angular. Velocity. 

( OSO.fSpacecr.af t ) ' - 

Ixy = *xz = ■’■yz 3 = 0.0; Ir xZ = 0.0 

. ± R „ - 0.05930; .0^(0) = 0. 03172 rad/sec. 
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Fig. 2k. Time History' of the Transverse Components of Main Body . Angular Velocity ' 

(SAS-A no damping) 

I xx = Izz = 27.00; Iyy = 2.8.54; Ixy = Iyz = 0; I xz = Y = 0.1048, w x (o) = 0.000159 rad/sec 
. = numerical integration; = approximate solution 
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Fig. :2 5 • Time History of the Transverse Components of Main Body Angular Velocity 
' (SAS-A No Damping) 

Z xx = I zz 3 2T.00; Iyy = 28. 5U; 1^ « I yz = 0; I xz = Y = 0.50j « x (0) = 0.0001 59' rad/sec 



TIME, t,- SECONDS 






!t ■ ^ - ** r ” + v 

y C&2 

2 + rtlf r 2 


r / ri J + ; 2 . 

_ 1 1 v* + A ? • 

^-s. 2 ■'*tf**i 

i CoSii r 

y 2 ~ C OS0 

Sin2 v * fr . + _ 3 ' Cos vj .♦ 


6r_ + 

° r J°o S 0 ,2 

V + rr 


~ i<U 

J C °S0 ;2 2 

3 + <r . 

° *i C °S0 ;2, 

^rr + V J + 

o ^ r. ~ i, ‘ 


1 T (r t * 

■2to u , ° rjC ° S ^) J * 2 * 

y w 2 4,f 2 + x 2 (<s- o 

* 1 2 + 2 W : , i (Si * 2 d, . 

y^l fr r p 1 9 3 + Sin 

r 2. ° r l c os0 + , 

V + 2,„ * 1 r .^ + o. • 


+ j-2, u * 

1 r r 1 2u> y*? (r r „ 

° 1 ° H 3 + r 2 ) ° l C °s<i> 
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4 
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2 


4-j 
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1 
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\ 

\ 


+ 
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la * aoj. 1 ! 19 ! 1 
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, 2 * r ! 9 2 v 1 
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g x 51 " ^ 1 * . , + i 4 si^4 

* A " <^3 siM3 
• - *2 COs4>2 . 
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B. Fortran Listing of Computer programming: 


// FORTRAN i '. . 

__*ONE WORD INTEGERS _ . _! = 

*LIST SOURCE PROGRAM ... 

SUBROUTINE SEN I (T,Y,DY) 

DIMENSION Y ( 6 ) » DY ( 6 > , C(5,5), CX{5,5) . . ..... 

COMMON WX, WY, WZ , S, Wl . PHI , A I XX, A I YY, A I ZZ , A I XY, AIXZ, AIYZ, AIRXX, 
$ AlRYYt AIRZZt'AIRXYt AIRXZ v AlRYZf AL>AK|RO*RLf AKK t AKf ALBAR t ALRY, 

$. RRY « AMBAR , N , M . . ; 

COMMON C,CPA,CPB, A16,A17, AL8, AL9, AllO, Alii, A112, 

SA113. AllAf AI15. AH6t A27.A28, A29,A2IOi A211, 

$A212,A213,A21A,A37,A38,A39,A310,A3U,A312,A313, 

$A31A,A315,A316,AA5,AA6,AA7,AA8, AA9 , AA l 0 , AA l 1 , 

$AA12,A5A, A55,A56,A57 

c ~ ' 7 7" ' 

DO 11 1=1, M 

DO 11 J = 1 , M __ _ _ . _ 

11 CX( I , J)=C( I , J) 

CX ( 1 , 2 ) = CX ( 1 , 2 ) -CPA*PHI 

CX ( 2 , 1 ) =CX(2, 1 )-CPA*PHI ! : . 

CX ( 3 , 1 ) =CX( 3, 1 )-CPB*PHI 

DY { 1 ) = (A16 + A17*PHI,)*WX*WY-A18*WY*WZ + A19*WY**2+( AUC*VsZ + 

........ $ A 1 1 1 * WY ) * S- A 1 1 2 W Z* * 2 - A 1 1 3* WZ^W X + A 1 1 A * W L **2 + A 1 1 5* >i 1 * « Y 

$All6*Wl*WZ 

DY ( 2 ) = ( A27+A28*PHI ) *WY* WZ-A29*WZ * WX-A2 1 0* WX*WY- 

$A21i*< WX**2-WZ**2)-( A212*WX-A213*-WZ ) * S + A 2 L 4 * PH l * W 1 * WY 

DY ( 3) = A3 7*V*X*WZ + A38*WX#*2-A39*WX4'WY- {A310+A311*PFI)=i=WY*WZ- 
$(A312+A313*PHI)*WY**2- ( A3 1 A*WY + A3 1 5*WX ) *S- A 3 1 6* PH I * W 1 * * Y 

DY ( A) = AA5-AA6*WZ*WX + AA7*WY*WZ + AA6-WZ*S + AA9*WZ* ! !'2-AA1C*WX**2- 

SAA11*WX*WY-AA12*WX*S 

DY ( 5 ) =A54*PHI 1 *'WY I ! tJ ! c 2 - A55 1 * : WX ! < t WY-A56 J>: WI-A57 ! <''PHI 
DY ( 6 ) = W 1 .. 

CALL SIMG(CX,CY ,M,KS) 

IF(KS) 3,2,3 

2. RETURN ... 

3 WR I T E ( 5 , A ) 

A FORMAT!//* SINGULAR EQUATIONS’) 

RETURN _ 

END ‘ 

FEATURES SUPPORTED _‘__1 

ONE WORD INTEGERS . • ** 

-CORE REQUIREMENTS FOR SEN1 .. 

COMMON 19A VARIABLES . 78 PROGRAM 512 



I 


PAGE' 2 V/28/72 “ 

RELATIVE ENTRY POINT ADDRESS IS 005 F (HEX) 

END OF COMPILATION 

\ T/~ "D U P . " ~~ • ■ . 

♦ STORE _ _ WS_ U A" SEN1 _ . 1 

CART ID 000 A ' DB AD DR 58 80 DB CNT 002 A 


// EJECT 


j 


“PAGE 3 4/28/72 " • “'1 

■ ■■■■■<-. -- - E 

* 

// FORTRAN _ _ _ • \ 

♦ONE WORD INTEGERS " ' ~ ?• 

♦LIST SOURCE PROGRAM , 

SUBROUTINE SEN 2 IT, Y, OY, IHLF, NDUMY » P) 

DIMENSION Y ( 6 ) * DY (6) , P (5) “ •' |j 

COMMON WX , WY , W Z , S » W 1 , PH I , A I XX , A I Y Y , A I Z Z , A I XY , A I X Z , A I YZ , A I R XX , P 
_ $_ AIRYY, AIRZZ, A IRXY, AI R X_Z_* A IR Y Z , AL , AM , RO , R1 t AKK , AK , AL BAR, A L RY, j 


$ RRY , AM8 AR , N » M 
DATA SY/O.O/ 
DATA WYL/O.O/ 


c 

c 

DATA TL/G.O/ 

■■ ■■ - I 

c 

SY= SY + O • 5* { T-TL ) * ( WY + WYL ) 
TL = T 



WYL =WY 

TOUT = T-S-0.C0C5 

CHECK =0.5* ( AIXX+AIRXX>*WX*WX+0. 

5 * A IYY^WY^WY+O. 5* 

$(AIZZ + A IRZZ )*WZ*WZ+ALRY*SY 
WRITE (5,4) TOUT, Y, IHLF, CHECK 
4 FORMAT ( 1X,F7.3,6E13.5, 13, E 13.5 ) 

P ( 5 ) = I N8 I T { 0 ) 

RETURN 

END 

FEATURES SUPPORTED 


[ ONE 

WORD INTEGERS 



CORE REQUIREMENTS FOR SEN2 
COMMON 58 VARIABLES 18 PROGRAM 130 


RELATIVE ENTRY POINT ADDRESS IS 001F (HEX) 

EN d“"o f “compTl a tTc n ~ 

// PUP . 

♦STORE WS UA SEN2 

CART I D O 0 OA DB ADO R 58 A A OB C NT O OOC 


// EJECT 
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PAGE 4 4 / 28/72 


/_/ FORTRAN _ _ i 

♦ONE WORD INTEGERS ? 

♦LIST SOURCE PROGRAM | 

♦I0CS(14C3 PRINTER) j 

EXTERNAL SENl » SEN 2 j 

DIMENSION C ( 5 » 5 ) , Y(6), DY(6), AUX (8,6), V(6) 1 

DIMENSION P ( 5 ) ! 

'COMMON' WX,WY,WZ,'S, WI , PHI ,AIXX,AIYY,AizZ,AIXY,AIX2,AIY'2,AIRXX, 

$ AIRYY,AIRZZ,AIRXY,AIRXZ,AIRYZ,AL,AM,RO,Rl ,AKK, AK , ALBAR , ALRY , - 


i 

i 




i 

♦ 

i 

i 


$ RRY,AMBAR,N,M ! 

COMMON C » CP A , CPB , A~l 6 , A 17 , A~1 8 , A 1 9 , A 1 10 , A 1 1 1 , A 1 1 2 , 

SA113,A114, A115, A116, A27,A28 , A29 , A2 1 0 , A2 1 1 , 
$A2I2,A213,A214,A37,A38,A39, A3 10 , A3 1 1 , A3 1 2 , A313, 

$ A 3 14, A 315, A3 1 6 , A 4 5 , A A 6 , A4 7 , A48', A 4 9", A 4 1 0 , A 4 1 1 , 

$A4 12 , A5A, A55, A56, A57 

EQUIVALENCE (P(l),TO), ( P ( 2 ) , TM ) , - ( P ( 3 ) , DT ) , (P(4),ERR) 

E Q U I V'A'L E NC E Ty ( 1) , W X ) , ( Y'f 2) , W'Y )', ( Y ( 3) , W Zl ,' ( Y ( 4 ) , ST ‘ 

EQUIVALENCE ( Y ( 5 ) , W 1 ) , (Y(6),PHI) 

DATA DT/0.1/ ■ 

* DA Y A 10/40 6.00/ - ‘ - [ 

DATA TM/420.00/ | 

DATA V/0. IE-3, 0.5,0. IE-3 , 0 . 2E3 , 0 . 3E-3 , 0 . 3E-2/ t 

DATA TOL/i • OE-4 / ' ' ~ " ;• 


[ C __ INITIALIZE 'COMMON^ 
5 C 

! WX = 0 . 94519 E -04 

i WY = 0 . 41851 E CO 


o n 


WZ = 0.24974E-03 
S = 0.19691E03 
PHI =0.0 
Wl =0.0 
AIXX =15.00 
A I YY = 28.00 
AIZZ ' =2.00. 

A I X Y = 0.0 

A I X Z = 0 . 0 

AI YZ '"= 0.0 
A I R XX = 5.575E-3 
AIRYY = 1I.519E-3 


AIRZZ 

= 5. 

575E-3 

AIRXY 

= 0. 

0 

AIRXZ 

= 0. 

0 

AIRYZ 

= 0. 

0 

AL = 

0.35 


RO = 

0.025 


R 1 =. 

203 


AM =0 

.0 


AK =0 

.0 


AKK = 

0.0 


AM B AR 

= 132 

.33 ■ 

ALBAR 

= 0. 

3499 

ALRY 

= 0.567E-2 

RRY = 

0.0 




COMPUTE COEFFEC I ENTS FOR SU8RCUT I ME ‘SEMI* 



f 
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^ f 


! PAGE 


A/28/72 


Cl 1,1 
C( 1,2 
CJ. 1,3 
C( i»A 
C< 1,5 
JC(2, 1 
C (2 i 2 
C<2, 3 
CJ_2,^ 
C ( 2,5 
C ( 3 » l 
CJ 3 t 2 
C (3 i 3 
C ( 3 » A 
C(3, 5 
C < A , 1 
C ( A » 2 
jC ( A i 3 
C ( A i A 
C ( A » 5 
C(5, 1 


C ( 5 » 2 
C ( 5 » 3 
CC 5 , A 
C ( 5 , 5 
CPA = 
CPD = 


= A I XX + A IRX'X + A.0*AM*AL**2 + 2.0*AM*{ R0 + R1 )**2 
= - ( AIXY+A IRXY ) 

= — ( A I X Z + AIRXZ) 


= -AIRXY 

= 0.0 

= - { A J X Y + AIRXYJ 

= A IYY +A IRYY +A.0*AM*(Rb 
= -( A1YZ+AIRYZ) 

= AIRYY 


+ R1 )**2 


) = AM*Rl*( R0+R1 ) 

) = ~(AIXZ+ AIRXZ) 

-( AIYZ+AIRYZJ 
A IZZ+AIRZ Z +A. 
-AIRYZ 

- ( AM*R 1*ALBAR ) 


) =_ 
) 

) = 

) = 


0*AM*AL**2 +2.0*AM* ( RO+R 1 )**2 


) = 
) = 
) = 
) = 
) = 
) = 
) 


= -AIRXY 
= AIRYY 
= -AIRYZ 
A I R Y Y 
0.0 
0.0 


AM*R1*<R0 + R1 ) 

= - ( AM*R1*ALBAR) 

= _0 . 0 ' __ _ 

" = AM*R 1 * *2 * ( l . 0 - AM/ A MB A R ) 
AM* AL*R 1 
AM*R1*(R0+R1 ) 


.1 - 
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/ 


A 310 
A 311 
A 312 
A 313 
A 3 l A 
A 315 
A 316 
A 45 = 


= A I X Z +AIRXZ 
=AM*R1*(R0+R1 ) 

-Mil +_ AIRZZ-AIYY-AIjm^j 

=' AIYZ + AIRY Z 

= AIRYY 

= A I R Y Z 

= AIYZ + A I RYZ 
= A I X Y +AIRXY 

=AM* AL*Rl _ , 

'= 2.0* AM * R 1 * A L BAR 
= 2.0 *Af^*Rl*( R0 + R1 ) 

: A I XY + A1RXY 
* A M * A L * R l ' 

=A I XX +AIRXX -Mil -AIRZZ 

=_ MY l + A I RYZ 

= A I X Z + A I RXZ 
= AIRYZ 

= AIRXY 

=" 2.0*Ar'*R0*Ri 
= ' AIYZ + AIRYZ 
: A I XY +AIRXY 

= A I V Y + A I R Y Y - A I X X - A I R X X 
= A I X Z + AIRXZ 

=AN*R1*(R0+R1 ) _ 

= A I X Y +AIRXY 
= AM*R1*ALBAR 

= A I RX Y 

=A I RYY • 

=2 . 0*AM*AL*R 1 

: ALRY +RRY 


4.0*AM*AL**2-2.0*AM*(R0+R1 )**2 | 


I 

r i 








i 


C THE FOLLOWING IS THE PROGRAM TO COMPARE 

C THE_ NUMERICAL INTEGRATION AND APPROXIMATE 

C SOLUTIONS (EQS. (135) £ (136)) 

C ARIME St BRIME ARE THE TRANSVERSE AND 

C POLAR MOMENTS O F INERT I A _OF _THE M A .1 N .. B C D Y....R E S P E C T IV ELY. 

C AK1 & AK 2 ARE THE CONSTANTS WHICH CAN BE RELATED 

C WITH THE INITIAL CONDITION 


f 

b 


I D. I .MEN S.I.Q N _.DA I J A 5.Lt 3) : 

• WR I TE ( 6 » 2 1 ) 

21 FORMAT! IX, T6, 'T' ,T15, »WX' ,T28, 'WZ' ) 

T T = 0 ; 4 • 

N D A T = 0 
9 T = TT 

„C.__ 

, f C LISTING OF DATA * * •. 

AK 1 = 0.000159 v . 

AK2 _= 2.03031 10 ; 

A = 1 A3 . 07850 
GAMMA = 0 . 10A 8 

A I RYY = 0.01 1 5_1 ? 

WY = 0.5 
P = -0.765A720 

0 1090 1 A 

C = 0.A916578 
ARIME = 27.005575 

BRJME_=. 28.551519 

AK 3 = ( GAMMA/AR IME ) *( WY+ ( P*'A I RYY ) / (0*BR IME ) ) 
ABSl = ((Q*C)/ARIME)/2.0 

ABS2 = DABS ( ABSl) . 

AROOT = 0. 5* ( 1/ABS2 ) **0. 5 
AKA = (. (AK3*AK1 )/,2.0)*AR00T 


T = TT 

All = -1. 0/(2. 0**0. 5) 

CALL SER 1 ( T, SUM 1 ) 

SSUM1 = SUM 1 
T = 0.0 

CALL_SER1{T,SUM1) 

All = AI 1*( (2.506-SSUM1 )-( 2 . 506-SUM 1 ) ) 

C 

c 

. — j : 

i 


V: 

h 
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l 


• 1 


T = TT 

A I 2 = -1. 0/(2. 0*0. 5) 

CALL SER2(T,SUM2J_ ; 

SSUM2 = SUM 2 
T = 0.0 

CALL _SER2( T, SUM2) 

A I 2 = A 1 2 * { ( -2 • 506-SSUK2 ) - ( -2 . 506-SUM2 ) ) 


XX = -1.0/2'. 0 
TA = ( TT+A ) **2 

JAA__=_A**2 £ _ _ ; 

A I 3 = -1.0/4. 0*'( (ABS2*TA)**XX-{ABS2*'tTaT**XX) 

C 

C 

F 1 = A K 4 * ( (DC0S( AK2 ) *AI l)-( DSIN( AK2) *AI 3) + (DSIN( AK2 )*AI 2 ) ) 

F2 = AK4*(-( DC0S( AK2 ) *A I 3 ) - ( DCQS ( AK2 ) *A I 2 ) + ( DS IN ( AK2 ) *A 1 1 ) ) 

F3=A K4*(-(0SIN( AK2 )*AI 1 ) - ( DCOS ( AK2 ) *A I 3 ) + ( DCOS ( AK2 ) * A I 2 ) )_ 

F 4 = A K 4 * ( ( DS IN ( AK2 ) * A I 3) + (DS IN ( AK2 ) *A 12) +"(OCOS ( AK2 ) * A I 1 ) ) 
WX=AK 1*DC0S ( (ABS1*TA)+AK2)+F1*DCOS(AOS1*TA)+F2*DSIN(ABS1*TA) 
__ WZ = AK 1 *D$ I N ( ( ABS1*TA )+AK2 )+F3*DC0S ( ABS 1 * T A ) +F4*DS I N { A BS 1 * T A ) 
WRITE (6,100) T T » W X , W Z 
100 FORMAT( lX t F7.3,E13.5,E13.5) 






NDAT = NDATt 1__ \ 

DAT ( ND A T , 1 ) = T f f 

DAT ( NDAT , 2 ) =WX 

DAT (NDAT, 3) =WZ 

TT = TT+1. 

I F ( TT-450. 5 ) 9,7,7 

_ 7 WR I TE ( 7 t 10 1 )_( ( DAT ( I * J )_, J 5 I ,_3 ).» I = 1 f NDAT ) 

101 FORMAT ( 20A4 ) - * 

CALL EXIT 

END 


h 


• % 


s 
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SUBROUTINE SERL(T,SUM1) 

IMPLICIT REAL*8 (A^H.CHZ ) 

SUM 1 = 0.0 
ARIME = 27.005575 

C = 0 .4916578 

0 = . 0 1 090 1 A . 

FF = (Q*C)/ARIME 

A = 143 . 078 5 0_ j 

YY =-1.0 

DO 10 M = 1,27"- 

FACT. .= ...1 . 0 ; 

MM = M-i 

IF (MM.EQ.O) GO TO 11 

AA_ = ( YY)** MM 

GO TO 4 
j 11 AA = 1 

4 A J = ( 2_*MMJ +0. 5 

A J J = M M + 0.5 
Y = FF* ( ( T + A ) **2 ) 


St 

5= 

Y**A J J 

Sll 

= 

Y * * M M 

S2 


AA*S1 

S3 

= 

S2/AJ 

MM 

s 

( 2*MM ) + 1 

DO 

1 

J = 1 , MM 

FAC 

T 

= F AC T* J 


1 CONTINUE 
S4=S3/FACT 

S5 = S4*S 1 1 

SUM2 = SUM2 + S 5 
10 CONTINUE 

RETURN 

END 










SUBROUTINE SER2(T,SUM2) 
IMPLICIT REAl*0( A-H.O-Z) 

SUM2__ = 0 .0 

C = 0.4916578 
ARIME = 27.005575 

Q = .01 0901 4 

FF = (Q*C)/ARIME 
A = 143.07850 

Y V = r 1.0 

DO 10 M = 1,27 
FACT = 1.0 

MM = M-i ___ 

I F ( MM. EO . 0 ) GO TO IT 
AA= (YY)**MM 

GOT 0 4 

11 AA = 1 

4 A J= ( 2*MM ) — • 5 

AJ 1 = MM- 0.5 

Y s FF*TlT+T)**2) 

51 = Y**AJ1 
S 1 1 = Y*#MM 

52 = A A * S 1 

' ‘ ? ■ _- r _ : ■. ^ ■■ S3 - S2/AJ 

IF (MM.EQ.O) MK = 1 

MM = MM*2.0 
DO 1 J = l, MM 

FACT = FACT*J 

1 CONTINUE 

S4=S3/F ACT 

S5 = S4*S11 

SUM 1 = SUM 1 + S 5 
10 CONTINUE 

" RETUR_N 

END 

//GO . SYS I N DD •* 

/* 


\ 


A 
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APPENDIX C : 

Nutation Damper Bias Angle due to the Offset of the Space- 
craft Mass Center in the Transverse Plane 

During the despin operation of the SAS-A spacecraft 

the nutation damper was observed to be biased by a small 

6 , 1 6 

amount off its expected equilibrium position. This 

bias angle was observed to diminish as the main body spin 
rate decreased. This phenomenon resulted from the actual 


offset of the spacecraft mass center in the transverse 
plane due to small errors in the final mass balancing 
prior to launch. 

This same bias in damper angle has been observed in 
the current numerical study (e.g. Figs. 5b, 6b, Jc) and can 
be explained by the fact that during the development of the 
equations of motion the lateral center of mass shift due to 
the damper motion was not included. In this appendix, the 
forces and moments acting on the damper mass in the trans- 
verse plane will be examined together with the bias angle 
from Fig. 5b and the center of mass offset displacement 
then calculated. 

Analy s i s 

(a) Offset of Center of Mass Due to Centrifugal Force: 

For static equilibrium of the damper pendulum, the 


C-l 


torque caused by the centrifugal force, is balanced by the 

restoring spring action. Referring to Fig. C-l and following 

Ref. 1 6 , the torque equation can be expressed: 

IN = mrio 2 sin a r, - Kd> n =0 (C.l) 

v 1 l 


where m = mass of damper pendulum 

K = torsion wire spring constant 

r = displacement of damper mass from actual mass 
center 

r^= damper pendulum length 
u)y= spacecraft spin rate 

(In Eq. (C-l), the effect of the Coriolis force has been 

neglected; it will be considered subsequently and shown to 

be a higher order effect for the SAS-A system). 

From Fig. C-l, sin a = d^_ 

r i ' 

so that Eq. (C.l) becomes: 

IN = mru 2 r 1 (_d ) - K <f> 1 = 0 (C.2) 

r i 

From the geometry of the figure. 


2 

r 

1 



d 


2 


r = r* + r cos g + f sin g 
1 o 

g = <j> 1 + a 


( C . 3 ) 


C-2 



and 


fcos3=d+r sin 3 ( C . *0 

o 

The exact solution for the spacecraft center of mass 
displacement f as a function ofr^, r^ s and to is compli- 
cated for large <j>^. It is noted in Eq. (C.2) that <J> 
varies directly as us 2 for constant r. The solution for f 

y 

in the case- of SAS— A can- -he- obtained — hy using small bias 
angles in the linearized equations. 

The following assumptions can thus be made: 
d , f << r ± 



a, 3, <(>2 are small and 3 * 

r = r n + r 
1 o 

Thus from Eqs. (C.2), (C.3) and (C.U) 


f = [ 


K 


m 


(rT+r ) w z 
x o y 


+ r 


o^l 


(C.5) 


Using the parameters given for the SAS-A satellite: 
m = 0.2158 kg 
K = 6.10 x 10" 5 nt-m 


rad 


r Q = 0.025 meter 
r^ = 0.203 meter 


C-3 


<J>^ =0.006 rad (obtained from Fig. 5b) 
a) = id ( 0 ) = 0.5 rad/sec 

J v 

and substituting these values into Eq. (C.5)» 
f = 1.52 millimeter 

For small angles Eq. ( C . 3 ) can be written: 

r = r^ + r Q cos <J>^ + f sin tj> -j_ ( C . 6 ) 

Substituting the values of r^ , r Q , f and <j>^ into Eq. (C.6) 
r = 0.228 meter 
From Eq. (C.4) 

d = f cos (f> 1 - r Q sin ^ 

After substituting the numerical values of f, r Q and <f>q the 
above equation yields, 

d = 0.0013701 meter 

The magnitude of the torque due to the centrifugal force about 
the damper hinge point is expressed: 

N = inroad 

c.f. y 

= 0.2158 x 0.228x ( . 5 ) 2 x 0.0013701 
= 1.68 x 10 ^ Newton-meter 
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(b) Calculation of the Coriolis Force and Torque : 

The velocity of a particle relative to a space fixed 

reference can be expressed as: 

V = V + wxr (C.T) 

s r 

where, V = velocity of the particle relative to a space 

S . 

fixed reference 

= velocity of the particle in a rotating frame 
with angular velocity u) relative to the fixed 
plane 

Again the acceleration of a particle relative to space can be 
obtained from the following equation. 


dt 


dV 


s 

dt 


+ wxV ( C . 8 ) 

b 


space rot 

Substituting the expression for V s from Eq. (C.T) into 
Eq. ( C . 8 ) 


dV, 


dt 


s pac e 


dt 


(V + wxr ) , + iox(V + mxr) . . ( C . 9 ) 

r rot r 


Using, V r = dr 

dt 


we obtain 


rot 


C-5 



dV 
s 

dt 


spac e 


= a + iiixr + 2ojxV + u)x(wxr) 
r r 


( C .10). 


where V is the velocity measured by an observer rotating 
r . 

with this system. The term | 2ajxY r | is the magnitude of the 
Coriolis acceleration. - 
From Fig. C-l 


r = (r Q + r x cos^Ji + (f + r 1 sinij^k (C.ll) 

Differentiating Eq. (C.ll) yields: 

= V r = ( -r-^s in^^i ) i + (f + r^ cos^^)]! . . ( C . 1 2 ) 

The total angular velocity vector during despin can be 
approximated by its largest component, toy as 

w = 'w y J (C.13) 


where 


a) is the angular velocity of the main body. 

«/ 

The Coriolis acceleration may be represented: 
*cor = 2[ “ xV r 3 


(C.lU) 


Substituting the values of w and into Eq . (C.lM 

a cor = 2 nuyj) x {(-r 1 sin<^ 1 4> 1 )i + (f + cos^ij )k } ] 



) 


and. after expansion, 

= [to (f + cos^^)]i - [ io (-r 1 sin<{> 1 4 ) 1 )]k (C.l6) 


cor 


The magnitude of the Coriolis force 


F cor = l F corl = 2 m/ “§C(f+ r^os^^)^ (-^ sin^^) 2 ' 


(C.17) 


For the static case f and 4>^ are both zero, so. 


F_ __ = 0 

cor 


Differentiating Eq. (C.5), for to = constant 

y 


^ m ( r + r )io£ + r o ^ ^1 
-L o «y / . „• * 


Substituting Eq. (C.l8) into Eq. (C.17) we obtain the 
result as , 


(C .18 ) 


/ 

F„ rtV< = 2m / a 2( 
cor (U . r 


m 


K . 

U + r, )( 0 Z + r o + r l cos ^) 2 ^ + r i sin 2 <|) 1 4.2 

O X y 


(c.19) 


From the time history of ^ for the case considered in Figs. 
5 an approximate average value of can be obtained as 


l*i 


ave 


= I 0.15 x 10 ^ rad/sec 


C-7 



into 


Substituting the SAS-A parameters and 4> 
Eq. (C.19), yields - 


1 av< 


F_ _ _ = 7.5^090 x IQ" 5 Newton 
cor 


The magnitude of the torque produced hy Coriolis 
force can be expressed by 


W cor = F cor cos a 


(C . 20) 


Using SAS-A nominal parameters it can be shown that: 

N = 1.53080 x 10 -6 Newton-meter, 
cor 


The average torque produced by Coriolis forces is about an 
order of magnitude less than the torque produced by the 
centrifugal force, so the effect of Coriolis force can be 
neglected in an approximate first order analysis. 


C-8 



Centrifugal Force Vector 
normal to pendulum 
= m r a) 2 sin « 

(balanced by torsion 
wire torque) 


Torsion Wire / 

Axis 

Pendulum 



FIGURE C-l: Geometry of Nutation Damper Bias Analysis 
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APPENDIX D 

Integrals Involved in the Solution of Equations (120) and 

( 121 ) 11 


■ y-1 

/ x sin ax dx 


ifCp-D 


2 a 4 


{ e 


F(y ,-iax) 


Il(i-y) 

2 

+ e r(p, iax) } Rey < 1 

a > 0 
x > 0 


fx 


P -1 

cos 


axdx 



y it 
2 


r (p, -iax) 


_iPTT 

+ e T(y , iax) } 



D-l 


